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ABSTRACT 

We present radial mass profiles within ~ 0.3 r v i r for 16 relaxed galaxy groups-poor clusters (kT 
range 1-3 keV) selected for optimal mass constraints from the Chandra and XMM data archives. 
After accounting for the mass of hot gas, the resulting mass profiles are described well by a two- 
component model consisting of dark matter (DM) , represented by an NFW model, and stars from the 
central galaxy. The stellar component is required only for 8 systems, for which reasonable stellar mass- 
to-light ratios (M/Lk) are obtained, assuming a Kroupa IMF. Modifying the NFW dark matter halo 
by adiabatic contraction does not improve the fit and yields systematically lower M/Lk- In contrast to 
previous results for massive clusters, we find that the NFW concentration parameter (c v i r ) for groups 
decreases with increasing M v j r and is inconsistent with no variation at the 3cr level. The normalization 
and slope of the c V i r -M V i r relation are consistent with the standard ACDM cosmological model with 
as = 0.9 (considering a 10% bias for early forming systems). The small intrinsic scatter measured 
about the c v i T -M v { r relation implies the groups represent preferentially relaxed, early forming systems. 
The mean gas fraction (/ = 0.05 ± 0.01) of the groups measured within an over-density A = 2500 
is lower than for hot, massive clusters, but the fractional scatter (it/// = 0.2) for groups is larger, 
implying a greater impact of feedback processes on groups, as expected. 

Subject headings: Cosmology: observations — dark matter — galaxies: halos — X-rays: galaxies: 
clusters — methods: data analysis 
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1. INTRODUCTION 

The properties of dark matter (DM) halos are a pow- 
erful discriminator between different cosmological scena- 
rios of structure formation. Dissipationless simulations 
of Cold Dark Matter (CDM) models find that the radial 
density profiles of DM halos are fairly well described be- 
tween approximately 0.01-1 r v j r (where r v ; r is the virial 
radius) by the 2-par ameter NFW model suggested by 
iNavarro et alj (|1997f ). Of particular importance is the 
distribution of halo concentration (c v i r , the ratio between 
r V j r and the characteristic radius of the density profile, 
r s ) and M v i r , the virial mass. Low mass halos are more 
concentrated because they collapse earlier than halos of 
larger mass, thus producing a predicted correlation be- 
tween Cvir and Mvi r . A significant scatter at fixed virial 
mass is expected and thought to be related to the distri- 
bution of halo forma tion epoch (e.g.. lBullock et aLll2001l : 
IWechsler et al.ll2002T ) . The c v - K -M v - lr relation and its scat- 
ter is a source of deviation from the self-similar scaling 
relation expected if the observable properties of halos 
are driven by simple gravitational collapse of the dom i- 
nant dark matter component (e.g.. iThomas et aI1l2001f ). 
For the currently favored ACDM model the median c v - lr 
varies slowly over a factor of 100 in M v j r , whereas the 
scatt e r remains very nearly constant (e.g., iBullock et al.l 
120011: iDolag et al.l 12004 IKuhlen et atll2005l ). The pre- 
cise relation between c v i r and M v j r is expected to vary 
significantly as a function of the cosmological parame- 
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ters, i ncluding g-g and w, the dark energy equation of 
state (IDolag et al.ll2004 IKuhlen et al.ll2005h . making an 
observational test of this relation a very powerful tool for 
cosmology. 

High quality X-ray data from Chandra and XMM 
observations indicate that the NFW model is a re- 
markably good description of the mass profiles of mas- 
sive galaxy cluste rs out to large porti o ns of their 
virial radii (e.g., iPratt fc Arnaudl l2002t iLewis et al.l 
20031: IPointecouteau et all 120051 : IVikhlinin et~ai] 120061 : 
Zapp acosta et al.l 12006ft . Typical values and scat- 



ter of concentrations determined from th e samp les of 
clusters analyzed in IPointecouteau et al.l (|2005h and 
IVikhlinin etail ([2006) are in general agreement with the 
simulation results. The Cy i r -M v j r relation measured by 
IPointecouteau et ail (|2005f K when fitted with a power 
law, has a slope of a = —0.04 ± 0.03. This slope is 
quite consistent with a constant value and is marg inally 
consistent (« 2a) with ACDM (IDolag et al.l l2004h. The 
optical study by lLokas et al.1 (|2006l ) using galaxy kine- 
matics for six nearby relaxed Abell clusters obtained re- 
sults consistent with the above X-ray studies but with 
larger uncertainty (e.g., a = — 0.6 ± 1.3). 

Observational tests of ACDM have proven controver- 
sial a t the galaxy scale (see discussion in lHumphrev et ahl 
|20 Q6D. Recently, usin g high quality X-ray Chandra data, 
in iHumphrev et al.l (|2006h we obtained accurate mass 
profiles for 7 elliptical galaxies, well described by a two- 
component model comprising an NFW DM halo and 
a stellar mass component. Omitting the latter com- 
ponent, which dominates the mass budget in the in- 
ner regio ns, leads to unphysica lly large concentrations 
(see also iMamon fc Lok as 2005) and may explain some 
large values found in the literature for elliptical galaxies 
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(jSato et all l200d : iKhosroshahi et ail |2004[ ). The mea- 
sured Cvir-Mvir relation of the 7 galaxies generally agrees 
with ACDM, provided the galaxies represent preferen- 
tially relaxed, earlier forming systems. 

Very few constraints exist on the group scale, where 
the simulations of DM halos are more reliable, com- 
pared to massive clusters, because a large number of ob - 
jects can b e simu lated at once fe.g- lBullock et al.ll2Q0l . 
iSato et all (|2000f ) investigated the c-M relation in X-rays 
using ASCA for a sample of objects ranging from 10 12 
to 10 15 M Q , including objects in the mass range dis- 
cussed in this paper. (However, neither the names of 
the objects in their sample, nor the description of the 
data reduction and analysis, has appeared in the litera- 
ture.) The slope obtained for the c-M relation was steep, 
-0.44±0.13. The limited spatial resolution of ASCA and 
energy dependence of its PSF made problematic the de- 
termination of reliable density and temperature profiles, 
and the authors neglected any stellar mass component in 
their fit s. Optical studies of grou ps using galaxy-galaxy 
lensin g (Mandelbaum et al.ll2006f) and caustics in redshift 
space iRines fc Diaferiol (|2006h obtain c v i r -Af v i r relations 
that are consistent with CDM simulations and with no 
variation in c with M, but with large errors. 

The scale of galaxy groups is also particularly intere- 
sting for the investigation of the influence of baryons on 
the DM profile. While the stellar mass component is 
clearly distinguished from the NFW DM component in 
the gravitating mass profiles obt ained from Chandra ob - 
servations of elliptical galaxies (Humphrey "eTaII[2r)M l. 
X-ray studies of relaxed clusters do not report signifi- 
cant deviations fr om a single NFW profile fitted to the 
gravitating mass dLewis et~a l. 2003; Pointccoutcau et all 
l2005t IVikhlinin et al.ll200llZappacosta et al.ll2006h. ex - 
cept for a few group-scale objects (Vikhli nin et al.ll2006f) . 
The group scale seems to represent a transition in the 
character of the mass p rofiles (and temperature profiles, 
iHumphrev et al.l I2006D and needs to be systematically 
explored. 

X-ray studies of mass profiles in galaxy systems have 
the advantage that the pressure tensor of the hot gas 
is isotropic and the gas in hydrostatic equilibrium (HE) 
traces the entire 3D cluster potential well. If one is care- 
ful to choose relaxed objects (i.e., with smooth, regular 
X-ray images) then hydrostatic equilibrium is a good ap- 
proximation and the resulting gravitating mass is reli- 
able, accurate to at least ~ 15% even in the presence of 
turbulence (e.g., lEvrard et al.lfl996t iFaltenbacher et all 
l2005t iNagai et al.ll2007f ). Because of limitations of previ- 
ous X-ray telescopes like ROSAT and ASCA, some sim- 
plifying assumptions like isothermality had t o be made 
for t he determination of group masses (see iMulchaevi 
l2000l and references therein) . Chandra and XMM have 
provided for the first time high quality, spatially re- 
solved spectra of the diffuse hot gas of X-ray groups (e.g. , 
Buote et al . 2003; Sun et al.ll20"oalO'Sullivan et al.| [2003l: 



Buote et al.ll2004l : iPratt fc Arnaudll2005l ). 



An investigation of the detailed mass profiles of galaxy 
groups (M = 10 13 -10 14 M ) with higher quality Chan- 
dra and XMM data is, therefore, timely. In this paper 
we present measurements of the mass profiles of a sam- 
ple of 16 groups chosen to provide the best mass de- 
terminations with current X-ray data. We selected the 
objects both to be the most relaxed systems (i.e., very 



regular X-ray image morphology), to insure hydrostatic 
equilibrium is a good approximation, and to have the 
highest quality Chandra and XMM data, which allow 
for the most precise measurements of the gas density 
and tem perature profiles. This paper is part of a series 
(see also IHumphrev et al1l2006t IZappacosta et al.l [20061 : 
iBuote et al.ll2006f) using high-quality Chandra and XMM 
data to investigate the mass profiles of galaxies, groups 
and clusters, placing constraints upon the c v i r -M v i r rela- 
tion over « 2.5 orders of magnitude in Af v j r . It is also 
the first in a series investigating the X-ray properties of 
groups and poor clusters: in future papers we will inves- 
tigate the entropy and heavy element abundance profiles. 
This paper is organized as follows. In §0 we discuss the 
target selection and in § [3] the data-reduction. We dis- 
cuss the spectral analysis in § [5j the mass analysis in § [6] 
and present the results in § [7] We discuss the results 
for individual objects in the sample in § [HI comparing 
with previous work in the literature. The systematic un- 
certainties in our analysis are discussed in § [9l and we 
present a discussion of our results in § [TU] with our con- 
clusions in § llll All distance-dependent quantities have 
been computed assuming H Q — 70 km s _1 Mpc -1 , fl m = 
0.3 and £1a= 0.7. Our assumed virial radius is defined as 
the radius of a sphere of mass M V j r , the mean density of 
which (for redshift 0) is 101 times the critical density of 
the universe (appropriate for the assumed cosmological 
model) and estimated at the redshift of the object. We 
will quote values for concentrations and masses at differ- 
ent over-densities to ease comparison with previous work 
in Appendix [XJ Our analysis procedure is described in 
greater detail in Appendix [B] All the errors quoted are 
at the 68% confidence limit. 

2. TARGET SELECTION 

For this study we choose, whenever possible, to focus 
on X-ray bright objects observed by both Chandra and 
XMM to exploit the complementary characteristics of the 
two satellites. The unprecedented spatial resolution of 
Chandra allows the temperature and density profiles to 
be resolved in the core, allowing us to disentangle the 
stellar and dark matter components. The unprecedented 
sensitivity of XMM ensures good S/N even in the faint 
outer regions, which is crucial because good constraints 
on the virial mass of the halo require density and temper- 
ature constraints over as large a radial range as possible. 
We looked for bright objects in the temperature range 
1-3 keV with sufficiently long exposures in the Chandra 
and XMM archives, together with our proprietary data. 
The potential targets were processed (§ [3]) and the im- 
ages in the 0.5-10 keV band examined (§[4| for evidence 
of disturbances: we choose objects which have a very reg- 
ular X-ray morphology, showing no or only weak signs of 
dynamical activity, with the peak of the emission coinci- 
dent with a luminous elliptical galaxy which is the most 
luminous group member. The only exception is RGH80 
which has two elliptical galaxies of comparable sizes in 
the core and probably a submerging group in the south 
(|Mahdavi et al.ll2~005( i. We include this object because 
it is part of a complete, X-ray flux-limited sample of 
15 groups that is scheduled to be observed by Chandra. 
It also allows an interesting comparison of derived mass 
properties with those obtained for the obviously relaxed 
systems in the sample. 



Dark Matter and Gas Fraction in Galaxy Groups. 



3 



TABLE 1 

The group sample and journal of observations 



Group 


z 


Dist 


AOlb 


Chandra exp 


T7DTP 


pn 


XMM exp 


c 

i out 








Mpc 


Aim point 


(ks) 


nltci 


mode 


(ksj 


kpc 




NOC ^044 


0090 


38.8 


S 


20.2 


Thin 


FF 


19 ^/IQ 3/8 9 4- 38 4/38 4/32 fl b 


326 


101.9 


NGC 1550 


0.0124 


53.6 


1 


9.8 _|_ g.6 a 


Medium 


FF 


21.4/22.6/17.8 


213 


102.2 


NGC 2563 


0.0149 


64.5 






A/Icdium 


FF 


20 4/20 8/16 5 


219 


102.4 


A 262 


0.0163 


70.7 


s 


28.7 


Thin 


EFF 


23.5/23.4/15.0 


254 


102.5 


NOC 533 


0.0185 


80.3 


s 


36.7 


Thin 


FF 


38 1 /37 4 /30 1 


271 


102.7 


MKW 4 


0.0200 


87.0 


s 


29.8 


Medium 


EFF 


14.0/13.9/9.4 


336 


103.1 


IC 1860 


0.0223 


97.1 






Thin 


FF 


34.1/34.8/28.0 


323 


103.1 


NGC 5129 


0.0230 


100.2 






Medium 


FF 


10.9/12.0/10.7 


241 


103.1 


NGC 4325 


0.0257 


112.2 


s 


30.0 


Thin 


FF 


20.8/20.8/14.7 


238 


103.3 


ESO 5520200 


0.0314 


137.7 






Thin 


EFF 


32.2/32.2/26.7 


418 


103.8 


AWM 4 


0.0317 


139.0 






Medium 


EFF 


17.5/17.2/12.5 


455 


103.9 


ESO 3060170 


0.0358 


157.5 


I 


13.8 + 13.9 a 






245 


104.2 


RGH 80 


0.0379 


167.0 


s 


38.5 


Thin 


EFF 


32.8/32.6/26.3 


533 


104.4 


MS 0116.3-0115 


0.0452 


200.2 


s 


39.0 






350 


105.0 


A 2717 


0.049 


217.7 






Thin 


FF 


49.2/49.6/42.9 


730 


105.3 


RXJ 1159.8+5531 


0.081 


368.0 


s 


75.0 






625 


108.0 



Note. — Listed above are the groups in our sample. Redshifts were obtained from NED and the distance is the inferred luminosity distance 
for a cosmological model with Hq— 70 km s~ 1 Mpc~ 1 ,Sl m — 0.3 and Ha— 0.7. The ACIS aim-point refers to S if the aim-point is located on 
the S3 chip or to I if the aim-point is located on one of the ACIS-I chips. The ACIS mode of all the observations was the Very Faint mode. 
The pn mode refers to FF if it is Full Frame or to EFF if it is Extended Full Frame; the MOS detectors were always in FF mode. The exposure 
times are net exposure times, after flare cleaning as described in the text, and for XMM they refer to MOSl/MOS2/pn. 
a observed twice with ACIS-I, with the core centered on ACIS-I0 in one occasion and on ACIS-I1 in the other. 
b the first set of exposures refers to the central pointing and the second set to the offset pointing. 
c the outer radius used in our analysis. 

d the assumed over-density, calculated at the redshift of the object. 



The details of the observations are given in Table [TJ 
We do not consider for analysis the available XMM ob- 
servations of ESO 3060170, MS 0116.3-0115 and RXJ 
1159.8+5531, because they are heavily contaminated by 
flares. We also do not consider the Chandra observa- 
tion of ESO 5520200 because of insufficient S/N for our 
purposes. In order to use as large a radial range as 
possible for objects observed in the ACIS-S configura- 
tion but lacking XMM data (MS 0116.3-0115 and RXJ 
1159.8+5531), we follow I Vikhlinin et alj (|2006l ) and also 
use the ACIS-S2 CCD in the analysis. 

The present sample has been selected preferentially for 
X-ray brightest and most relaxed groups to obtain the 
best constraints on the mass profiles in individual ob- 
jects with current data. Although the sample is biased 
and is not statically complete, our analysis of these sys- 
tems represents an essential step in the investigation of 
DM in galaxy groups with X-rays. In future work we 
will compare these results to those obtained using the 
complete, X-ray flux-limited sample of 15 groups noted 
above. 

3. DATA REDUCTION 
3.1. Chandra 

For data reduction we used the CI AO 3. 2 and Hea- 
so/i 5.3 software suites, in conjunction with the Chan- 
dra Caldb calibration database 3.0.0. In order to en- 
sure the most up-to-date calibration, all data were re- 
processed from the "level 1" events files, following the 
standard Chandra data-reduction threads . We applied 
corrections to take account of a time-dependent drift in 
the detector gain and, for ACIS-I observations, the ef- 
fects of "charge transfer inefficiency" , as implemented in 
the standard CI AO tools. From regions of least source 
contamination of the CCDs we extracted a light-curve 

4 http:/ /cxc. harvard.edu/ciao/threads/index. html 



(5.0-10.0 keV) to identify periods of high background. 
Point source det ection was performe d using the CI AO 
tool wavdetect ([Freeman et al.ll2002l ). The source lists 
created in different energy bands (so as to identify un- 
usual soft or hard sources) were combined, and dupli- 
cated sources removed. The final list was checked by vi- 
sual inspection of the images. The resolved point sources 
were finally removed so as not to contaminate the diffuse 
emission. Further de tails about the Chan d ra dat a reduc- 
tion can be found in iHumphrev k. B uotc ( 2006) . 

3.2. XMM 

We generated calibrated event files with SAS v6.0.0 
using the tasks emchain and epchain. We considered 
only event patterns 0-12 for MOS and for pn. Bright 
pixels and hot columns were removed by applying the 
expression (FLAG == 0) to the extraction of spectra 
and images. We correct statistically for the pn out-of- 
time (OoT) events. Following the standard procedure, 
we generate an OoT event list, processed in the same way 
as the observation, and then subtract it from the images 
and spectra, after being multiplied by the mode depen- 
dent ratio of integration and read-out time (6.3% for Full 
Frame and 2.3% for Extended Full Frame). The energy 
scale of the pn over the whole spectral bandpass has been 
further improved using the task epreject. We clean the 
data for soft proton flares using a threshold cut method 
by means of a Gaussian fit to the peak of the histogram 
of the 100s time bins of the light curve (see Append ix A 
of lPratt fc Arnaudl 120021 : iDe Luca k Molendilbooi ) and 
excluding periods where the count rate lies more than 
3cr away from the mean. The lightcurves were extracted 
from regions of least source contamination (excising the 
bright group core in the central 5' and the point source 
list from the SOC pipeline, after visual inspection) in 
two different energy bands: a hard band, 10-12 keV for 
MOS and 10-13 keV for pn, and a wider band, 0.5-10 
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Fig. 1. — X-ray images of each of the objects in the sample (see £|4t . 

keV, as a safety check for p ossible flares with soft spec- 
tra (Hiyiiiinmiran [2001 Tl^^A^m^^QQ^) ■ The 
flaring periods thus determined where further checked 
by visual inspection of the light curves. Point sources 
were detected using the task ewavelet in the energy band 
0.5-10 keV and checked by eye on images generated for 
each detector. Detected point sources from all detec- 
tors were merged, and the events in the corresponding 
regions were removed from the event list, using circular 



regions of 25"radius centered at the source position. The 
area lost due to point source exclusion, CCD gaps and 
bad pixels was calculated using a mask image. Redistri- 
bution matrix files (RMFs) and ancillary response files 
(ARFs) were generated using the SAS tasks rmfgen and 
arfgen, the latter in extended source mode. Appropri- 
ate flux-weighting was performed for RMFs, using our 
own dedicated software, and for ARFs, using exposure- 
corrected images of the source as detector maps (with 
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pixel size of 1', the minimum scale modeled by arfgen) 
to sample t he variation in emi s sion, following the pre- 
scription of ISaxton fc Siddiquil (|2002D . Spectral results 
obtained using ARFs are completely consi stent with the 
other frequently employed method (e.g. lArnaud et al.l 
1200 lh of correcting directly the spectra for v ignetting 
(|Gastaldello et alj|2003t iMorris fe Fabianll2005D . 

3.3. Background subtraction 

Insuring proper background subtraction is one of the 
key challenges associated with the spectral fitting of low 
surface brightness, diffuse, X-ray emission. The back- 
ground experienced by both Chandra and XMM consists 
of (1) an extreme time variable component due to soft (E 
~ tens of keV) protons channeled by the telescopes mir- 
rors, (2) a slowly changing (with variability time scale 
much longer than the length of a typical observation) 
quiescent component due to high energy particles (E > 
a few MeV), and (3) the sky X-ray background, decom- 
posed into the extragalactic Cosmic X-ray ba ckground by 
AGN and the Galactic X-ray emission (e.g. iLumb et al.l 
2002; IMarkevitch e"t~aT][2003fi . 

The "blank fields" distributed by the observatories are 
not a perfect representation of the background in any one 
observation. Firstly, there are significant long term vari- 
ations in the quiescent particle background. Secondly, 
the soft Galactic background component varies strongly 
from field to field. Finally, there may be some residual 
mild flaring. 

Two approaches have been investigated to obtain 
more accurate background estimates than provided by 
the "blank field" background templates: the double 
subtraction techniq ue (see details in Appendix A of 
lArnaud et alj I2002D and a complete modeling of the 
various background co mponents fe.g. ILumb et al.l [20021 : 
IMarkevitch et all [2003). Double subtraction is, in prin- 
ciple, very effective, but particular care has to be taken 
to locate a region in the field of view of the observa- 
tion completely free of source emission; this is difficult 
for nearby objects. The complete modeling of the vari- 
ous background components can rely on a large number 
of observations performed to characterize the quiescent 
particle background component (stowed or Dark Moon 
for Chandra, closed for XMM) , and on the large number 
of observations which constitute the "blank field" data 
sets to characterize the sky background components. The 
drawback is that the resulting model, which also includes 
a source component, is complicated, and parameter de- 
generacies can arise. However, the method is particu- 
larly effective for studying groups and poor clusters be- 
cause the source component (mainly characterized by the 
~ 1 keV Fe-L shell blend) is clearly spectrally separated 
from all the other background components. For the use 
of this approach to the Chandr a data we refer the reader 
to iHumphrev fc Buotd ([20061) . Here we will describe 
the procedure used for XMM data which elaborates an d 
updates the procedure described in iBuote et alj (2004). 
The algorithm implemented has the following main steps: 

- Characterization of the quiescent particle back- 
ground. We co-add individual spectra taken from 
closed observations. The spectra in the 0.4-12 keV 
band for MOS and 0.4-13 keV band for pn can be 
adequately described by a broken power law con- 



tinuum and several Gaussians for the instrumen- 
tal lines. Typical values for the model parameters 
are: 0.7-0.8 for the slope at low energies, 1.0-1.2 
keV for the break energies, 0.2 for the high energy 
slope for MOS, 0.4-0.5 for the high energy slope for 
pn. While the low-energy slope exhibits significant 
variation between the observations in our study, the 
high-energy slope is very stable. These results are 
consistent with previous s t udies (ILumb et alj|2002t 
iDe Luca fc Molendil 120041 iNevalainen et al.ll2005f) . 
The spectral shape of the continuum broken power 
law does not change significantly across the detec- 
to r, nor does it vary in time (as in the MOS study 
bv IDe Luca fc Molendill200l at high energies. 

- The model derived from the closed data is fitted 
to the spectra of the Out of Field of View (OFV) 
events of each observation. Portions of the MOS 
and pn detectors are not exposed to the sky, and 
therefore neither cosmic X-ray photons nor low en- 
ergy particle-induced events (like from soft pro- 
tons) are collected. Indeed, while this is almost ex- 
actly true for the MOS (the fraction of OoT events 
is 0.35% in FF mode), for the pn a higher frac- 
tion of in FOV events events are assigned to the 
OFV region as OoT (6.3% in FF mode and 2.3% 
in EFF). In the case of strong flaring, this OoT 
contribution can seriously affect the pn OFV spec- 
trum. We therefore chose to extract OFV events 
for the pn after the flare cleaning. The model de- 
rived from the OFV data is taken as the initial rep- 
resentation of the quiescent particle background for 
the particular observation. 

- To the broken power-law plus Gaussian lines de- 
scribing the quiescent instrumental background 
(not vignetted and implemented as a background 
model in Xspec) , we add the componen ts describing 
the sk y X-ray background, following ILumb et al.l 
(2002): a power law with s lope T = 1.41 and 
norma lization as reported in IDe Luca fc Molendil 
|[200j), free to vary within the cosmic variance as 

(|Barcons et al.l2000D where fl is the solid an- 
gle covered by the observation; two thermal compo- 
nents with temperatures 0.07 and 0.20 keV respec- 
tively, and abundances fixed at solar. When model- 
ing sources projected toward the North Polar Spur 
(NGC 5129) we found it necessary to add a third 
ther mal component at ~ 0. 4 ke V, in agreement 
with IMarkevitch et all (|2003l ) and iVikhlinin et all 
(|2005l) . 

- This model, plus a source component described by 
a thermal plasma with temperature and abundance 
free to vary, was fitted jointly to the outermost 
annuli used in the spectral extraction (see below). 
The parameters of the source component were free 
to vary in each annulus. The normalizations of 
the cosmic components and of the broken power 
law component were scaled according to the ratio 
of the annuli area, while the normalizations of the 
instrumental lines were free to vary, given the fact 
that these components are highly spatially variabl e 
(e.g. IDe Luca fc Molendil l200l ILumb et al J 120021) . 
An inter-calibration constant free to vary between 
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0.9 and 1.1 was added to the model to take into ac- 
count any cross-calibration differences between the 
three EPIC instruments. Given the best fit model 
for these annuli, we generate a pulse height ampli- 
tude (PHA) correction file used in Xspec. 

- For the annuli not involved in the background fit- 
ting, we scale the resulting model to the area of 
the annular region of interest in the spectral ex- 
traction and generate a PHA file. To take into 
account the variable instrumental lines, we renor- 
malized the instrumental line components in the 
model using the corresponding regions extracted 
from the background templates. 

We mention that possible slight variations in the par- 
ticle continuum across the dete ctor plane (see Appendix 
A of lDe Luca k Molendll2004D . or residual mild flaring, 
has been modeled by slightly changing the high energy 
slope of the broken power law. This does not have any 
tangible effect on the spectral parameters derived for soft 
X-ray sources like the objects considered in this paper. 

4. X-RAY IMAGES 

The X-ray image of each group was examined to iden- 
tify any significant surface brightness disturbances in- 
dicating departures from hydrostatic equilibrium. Low 
level X-ray disturbances like the weak signs of AGN ac- 
tivity in the center of Abell 262 (jBlanton et al.ll2004D . or 
the presence of a submerg ing group in the south of RGH 
80 (jMahdavi et al.l l2005[ h do not seriously impact X- 
ray mass determinations, provided care is taken to avoid 
highly disturbed emission (|Buote k Tsailfll9^[Schlndien 
1996). The images for 15 objects in the sample are shown 
in Fig[TJ for objects which have XMM data we show 
the combined MOS1 and MOS2 image in the 0.5-2.0 keV 
band. For those objects with only Chandra data (3 out of 
16) we did the following: For RXJ 1159.8+5531 and MS 
0116.3-0115 we display the 0.5-10 keV ACIS-S3 image, 
while for ESO 3060170 we show the 0.5-10 keV ACIS- 
I image. The images were processed to remove point 
sources using the CI AO tool dmfilth, which replaces pho- 
tons in the vicinity of each point source with a locally es- 
timated background. The images where then flat-fielded 
with a 1.7 keV exposure map for Chandra images and 
a 1.25 keV exposure map for XMM. Then we smoothed 
the images with a 5" Gaussian for Chandra and a 16" 
Gaussian for XMM to make large-scale structure more 
apparent. For both the Chandr a and XMM image s of 
NGC 5044 we refer the reader to lBuote et ail (|2003[ ). 

None of the objects show obvious large scale distur- 
bance in their X-ray emission. The only notable sub- 
structure is the infalling subgroup in the southern region 
of RGH 80, evident as a tail of enhanced emission. Some 
disturbance is also present in the core of RGH 80 as re- 
vealed by our Chandra image. We masked in our analysis 
the region of enhanced X-ray emission associated with 
the subgroup. For other systems some low-amplitude, 
small-scale disturbances are present, such as a surface 
brightness discontinuity, reminiscent of a cold front, in 
the NW of IC 1860; the cavities in the central 10 kpc of 
Abell 262, as reveale d by the Chandra image presented in 
iBlanton et alj (|2004h ; and the fi lamentary structur es and 
possible cavities in NGC 5044 (jBuote et all 12003). The 



"cooling wake" discussed i n the XMM image of NGC 
5044 bv iFinoguenov et alJ (|2006h . is simply the bright 
SE arm of the finger-like structure caused by the cavi- 
ties. Hint of caviti es have been detected in NGC 4325 
([Russell et al.ll2007h and there are signs of sloshing in the 
core of MKW 4. We asses the impact of these features 
in §iH 
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Fig. 2. — Temperature profile for the NGC 533 group derived 
from XMM data (black) and from Chandra data (red) . The inner 
two XMM bins have not been considered in the derivation of the 
mass profile. 



5. SPECTRAL ANALYSIS 

We extracted spectra in concentric circular annuli lo- 
cated at the X-ray centroid computed within a radius 
of 30", with the initial center on the peak of the X-ray 
emission. The widths of the annuli were chosen to have 
approximately equal background-subtracted counts and 
to have a minimum width of 60" for XMM to avoid under- 
sampling of the PSF. For Chandra, given the better PSF, 
the widths of the annuli, in practice, were only limited 
by count statistics. The spectra were re-binned to ensure 
a signal-to-noise ratio of at least 3 and a minimum 20 
counts per bin (necessary for the validity of the \ 2 mini- 
mization method) . We fitted the background-subtracted 
spectrum with an AP EC thermal plasma mo dified by 
Galactic absorption (iDickev k Lockmanl I1990D to each 
annulus. The free parameters are temperature, nor- 
malization (proportional to emission measure), and the 
abundances Fe, O and, when possible, Si and S. The im- 
pact of unresolved point sources, in particular LMXB in 
the central galaxy, was taken into account by adding a 
7.3 keV bremsstrahlung component for all annuli within 
the twenty-fifth magnitude isophote (D25) of the cen- 
tral galaxy, taken from LED A. (This model gives a good 
fit to the composit e spectrum of th e detected sources 
in nearby galaxies, llrwin et al J 120031 ) This is particu- 
larly relevant for the inner XMM annuli, where in general 
point sources are not detected. T he spectral fit ting was 
performed with Xspec (ver 11.3.1. lArnaudl[T996l ). We es- 
timated the statistical errors on the fitted parameters by 
simulating spectra for each annulus using the best fitting 
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TABLE 2 

Optical properties of the group central galaxy 



Name 


Group 


r e (B) r e (K) L B L K 






kpc (arcsec) kpc (arcsec) 10 10 Lq 10 11 Lq 



NGC 5044 


NGC 5044 


NGC 1550 


NGC 1550 


NGC 2563 


NGC 2563 


NGC 708 


A 262 


NGC 533 


NGC 533 


NGC 4073 


MKW 4 


IC 1860 


IC 1860 


NGC 5129 


NGC 5129 


NGC 4325 


NGC 4325 


ESO 552-020 


ESO5520200 


NGC 6051 


AWM 4 


ESO 306-017 


ESO3060170 


MCG 6-29-77 


RGH80 


MCG 6-29-78 


RGH80 


UGC 842 


MS 0116.3-0115 


ESO 349-22 


A 2717 


2MASSX J11595215 


RXJ 1159.8+5531 





4.53 (24.5) 


6.98 


2.87 


6.45 (25.5) 


3.05 (12.1) 


4.33 


2.09 


5.89 (19.3) 


3.73 (12.2) 


3.84 


2.66 


25.60 (77.1) 


10.16 (30.6) 


3.84 


4.12 


16.92 (45.4) 


9.22 (25.2) 


12.4 


6.14 


19.24 (47.5) 


10.25 (25.3) 


13.7 


7.18 


8.34 (18.5) 


8.03 (17.8) 


6.08 


4.38 


13.34 (28.7) 


6.60 (14.2) 


12.0 


4.99 




5.22 (10.1) 


4.61 


2.33 




15.7 (25.0) 


15.6 


8.19 


10.21 (16.1) 


10.33 (16.3) 


9.91 


7.50 




18.51 (26.0) 


18.5 


6.95 




5.11 (6.8) 




2.93 








6.23 (8.3) 


4.21 


2.39 




9.69 (10.9) 


9.29 


5.77 




15.53 (16.2) 


9.20 


5.42 




9.77 (6.4) 


23.6 


10.3 



antral 

r e in the B band from RC3 (dc Vaucouleurs ct al. 1991), Lk (using the luminosity distance of Table m 
and r e in the K band from 2MASS. 



models and then fitted the simulated spectra in exactly 
the same manner as done for the actual data. From 20 
Monte Carlo simulations we compute the standard devi- 
ation for each free parameter, which we quote as the "1 
ct" error (these error estimates generally agree very well 
with those obtained using the standar d A\ 2 approach in 
Xspec, e.g.. iHumphrev fc Buotell2006h . 

If an object has been observed by both Chandra and 
XMM we selected for our final analysis only the Chandra 
data in the inner core region were the temperature rises 
outward from the center. The XMM spectra extracted in 
wide annuli are not as well fitted by a single temperature 
emission model as are the Chandra spectra in narrower 
annuli, suggesting that departures from single temper- 
ature emission in the projected spectra stem primarily 
from the steep radial temperature gradient present in 
the core, as shown in Fig[2] for NGC 533. The better 
Chandra PSF also allows us to exclude point sources un- 
detected with XMM, in particular LMXB in the central 
galaxy. Unresolved LMXB still affect the Chandra data, 
but to a much lesser extent than XMM data. This com- 
ponent is evident as an excess at energies greater than 
~ 3 keV and can contribute, if neglected, to the multi- 
phase appearance of XMM spectra. 

6. MASS MODELING 

To calculate the gravitating mass distribution we solve 
the equation of hydrostatic equilibrium assuming spher- 
ical symmetry. By requiring spherical symmetry we ob- 
tain spherically averaged mass profiles which allows us to 
test the spherically averaged DM profiles obtained from 
cosmological simulations and to facilitate comparison to 
previous observational studies. 

Following the approach adopted in IHumphrev et all 
(2006), we assume parametrizations for the temperature 
and mass profiles to calculate the gas density assuming 



Mdr 
r 2 T 



(1) 



hydrostatic equilibrium, 

T ( -pm A G 

where r is the radius from the center of the gravitational 
potential, p g is the gas density, p g o and To are density 
and temperature at some "reference" radius ro, k B is 
Boltzmann's constant, G is the universal gravitational 
constant, tua is the atomic mass unit and p (taken to be 
0.62) is the mean atomic weight of the gas. The p g (r) 
and T(r) profiles are fitted simultaneously to the data to 
constrain the parameters of the temperature and mass 
models. Since the gravitating mass also contains the gas 
mass, eq. ([T]) is solved iteratively for p g . 

This "parametric mass method" is the principal ap- 
proach employed in this study. We assess systematic 
errors in this adopted method in ^9.51 by comparing to 
results obtained using other solutions to the hydrostatic 
equilibrium equation. Firstly, rather than solving for the 
gas density, we can solve for the temperature, 



T(r) = T 



Pgo 

p g (r) k B p g {r) 



piiriAG f r p g Mdr 



(2) 



which provides an alternative implementation of the 
"parametric mass method" . Note that in both cases the 
parameters of the mass model are obtained from fitting 
the gas density and temperature data. The goodness- 
of-fit for any mass model (e.g., NFW) can be assessed 
directly from the residuals of the fit. Secondly, we use 
the more traditio nal formulation of the hydrostatic equi- 
librium equation (Mathcwsl ll978f ). 



AI(< r)=r 



k B T(r) ( din 
Gpmj, 



dlnT 



d\nr d\nr 



(3) 



which involves parametrizing independently p g and T 
using simple functional forms in order to evaluate the 
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derivatives in Eq. ([3]). Since, however, the mass pro- 
file itself is not parametrized, we denote this traditional 
approach a "non-parametric mass method" . Since the 
mass profile itself is produced by this method, if one 
wants to evaluate the success of a particular mass model 
(e.g., NFW) then additional fitting is req uired. Con- 
seque ntly, following previous studies (e.g., iLewis etlrtl 
|2003|) for each annulus we assign a three dimensional ra- 
dius value r = [(R%1 + i?f n /2 )/2] 2 / 3 , where R in and R out 
are respectively the inner and outer radii of the (two- 
dimensional) annulus. At each radius r we calculate the 
total enclosed gravitating mass M(< r) according to the 
equation of hydrostatic equilibrium. The errors on the 
resulting mass "data points" were estimated from the 
Monte Carlo simulations used to estimate the errors for 
density and temperatures (SjSJ, giving a set of mass val- 
ues at each radius. From those we calculate the stan- 
dard deviation which we quote as the "lcr" error for this 
method. To analyze the shape of the mass profiles, we 
fitted parametrized models to the mass values. 

There are several reasons why we adopt equation [Tj 
instead of equation [3] for our analysis. Firstly, as noted 
above, the "parametric mass method" allows a partic- 
ular mass model to be constrained immediately by the 
gas density and temperature data and the goodness-of- 
fit of the mass model can be assessed in a straightfor- 
ward manner. Secondly, despite the high-quality X-ray 
data provided by Chandra and IMM, it is still not pos- 
sible to compute accurate derivatives of the tempera- 
ture and density profiles at each radius. Consequently, 
smooth models must be fitted to the entire radial profiles, 
which may not produce physical solutions to equation [3] ; 
e.g., jagged, non-monotonically increasing mass profiles. 
Thirdly, we analyze the projected temperatures and den- 
sities which requires the models for the gas density (den- 
sity squared, see below) and temperature to be projected 
along the line-of-sight. This requires evaluating the mod- 
els at least out to the virial radius, well outside the outer 
radius of most of the groups in our sample. By using 
the "parametric mass method" any extrapolation of the 
gas density (equation [1]) or temperature (equation [2]) 
is performed consistently within the context of the as- 
sumed mass profile. It is for this last reason we have 
a slight preference for using equation Q] over equation 
[2] for this study. Nevertheless, despite these differences, 
we find that the different approaches to the mass mod- 
eling represented by the three equations give consistent 
results, within the errors, for the global halo properties 
(see also ^9.5[) . 

For our default analysis we projected parametrized 
models of the three-dimensional quantities, p 2 g and T, 
and fitted these projected models to the results obtained 
from our analysis of the data projected on the sky (see 
$5]). The models have been integrated over each radial 
bin (rather than only evaluating at a single point within 
the bin) to provide a consistent comparison. They also 
have been projected along the line of sight including the 
radial variation in the plasma emissivity A(T, Zp e ), using 
a model fitted to the observed Zp e profile. We provide a 
review of the projection of spherical coronal gas models 
for comparison to X-ray spectral data in Appendix [B] 

On the left panels of FigEJ FiglH FigE] and FigE] we 
show the radial profiles of the emission- weighted projec- 



tion of p g 2 (i.e., proportional to the norm parameter of 
the APEC model divided by the area of the annulus) 
along with the best-fitting model and residuals; in the 
central panels we show the radial profiles of the measured 
T along with the best fitting emission- weighted projected 
model and residuals. 

6.1. Gas density models 

We considered two model s for fitting of the gas den- 
sity profile: the (3 model (Cavalie re fc Fusco-Femianol 
1978ft and a cus ped (3 model ([Pratt k. Arnaudl l2002t 
Lewis et all 12003ft . the latter of which is a modified 
model allowing for steepening of the profile in the inner 
regions (r < r c ). This model was introduced to account 
for the sharply peaked surface brightness in the centers of 
relaxed X-ray systems. This mo del has now been widely 
used for both low-redsh ift (e.g. iPratt fe; Arnaud 12002ft 
and high-redshift (e.g., iKotov fc Vikhlininl 12005 ) clus- 
ters. It is preferred wi th respect to the double-/? model 
(e.g.. iMohr et al.lll999f ) because the two models give fits 
of comparable quality, while the cusped (3 model has one 
less free parameter. The cusped (3 model is also bet- 
ter behaved in the mass determination using the "non- 
parametric method" defined above (eqn. [3| . 

6.2. Temperature models 

The projected temperature profiles for our groups 
show a large degree of similarity. We adopted several 
parametrizations that have enough flexibility for each 
system to describe the temperature profile reasonably 
well and to explore the sensitivity of our results to the 
particular functional form. The analytic models we con- 
struct are the following: 

- Smoothly joined power laws: 



T{r) = ■ 



1 



+ (t2(r)) 



i = 1,2 



lOOkpc, 

- Power laws mediated by an exponential: 
T(r) = T + ti(r)e _( ^ r + t 2 (r) ( 1 - e~ ( ^ } 



(4) 



U(r) = Ti 



i = 1,2 



(5) 



The lAllen et all (f200lh rising profile joined to a 
falling temperature profile by an exponential cut- 
off, 



T{r) = h{r)e ( "* V + t 2 (r) ( 1 - e 



h(r)=a + T 1 



t 2 (r) = b + T 2 



in 







1 + 













i + 


(ft)" 
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The third ( "RiseFall" ) model has been adopted in par- 
ticular for temperature profiles showing an inner core 
flattening like NGC 533, NGC 4325 and NGC 5044, while 
the first two models provide comparable fits to the gen- 
eral profile. We will assess how different choices of tem- 
perature profile, together with density profiles, affect our 
mass measurements in H9.5I 

6.3. Mass models 

We compute the total gravitating mass as the sum of 
DM, stars, and hot gas: Mdm + M star s + M gas . For this 
study we only consider the contribution of the central 
galaxy to the stellar mass. The X-ray data provide a 
direct measurement of the hot gas density and therefore 
of M gas . We tested the following mass models against 
the data: 



all of the mass components.) The values thus obtained 
have then been converted to various over-densities (in 
particular the virial over-density, listed for each object 
in Table [H) in Table UJ b y using the formula provided 
bv lHu fc Kravtsovi (|2003f ) appropriate for an NFW halo. 
We prefer this procedure for extrapolating the mass and 
concentration to A sa 101 (for comparison with theoret- 
ical models) because it does not involve also extrapolat- 
ing M gas . We find that the extrapolated values for the 
gas mass are sensitive to the radial range over which the 
density profile is fitted (see ^9.7|) . 



- M DM = NFW, _ M stars = 0: A single NFW 
model to inves ti gate scena rios like the ones o f 
ILoeb fe Peeblei (f200l and lEl-Zant et all (120041 ) 
and the effect of the omission of the stellar mass, if 
present, on the derived concentration parameter. 

- M DM = NFW, M stars = deV. A NFW model plus 
a model for a stellar component. We adopted a 
de Vaucolueurs stellar mass potential using the 
analytical approximation to the dep rojected Sersic 
model of Prugnie l fe Simie"nl ([1997) with n = 4. 
The de Vaucolueurs profile is a good description of 
the stellar light distribution even for objects which 
follow the more general Sersic profile with Sersic in- 
dex n ^ 4 (see appendix A of iPadmanabhan et al.l 
2004). The de Vaucolueurs effective radius r e is 
measured in the K-band by the Two Micron All 
Sky Survey (2 MASS ) as listed in the Extended 
Source Catalog (jJarrett et al.ll2000L see Tab©. We 
refer to this model as NFW+stars. 

- M DM = NFW* AC, M stars = deV. A NFW compo- 
n ent modified by the adiabatic contraction model 
of lGnedin et alj (|2004h 5 plus a de Vaucolueurs com- 
ponent for the stellar mass, to explore the impor- 
tance of baryon condensation in the central galaxy 
for the DM halo profile. We refer to this model as 
NFW*AC+stars. 

- Finally we also ex amined the recently suggested 
Sersic-like profile (|Navarro et al.l 120041 hereafter 
N04) which should be a better parametrization of 
the innermost regions of CDM halos. 

To obtain the true virial radius and virial mass (and 
concentration) we initially take r& and Ma obtained for 
the DM component, where A corresponds to the over- 
density level (2500, 1250, 500) closest to the radial range 
covered by the data and listed in Tabled! Then we added 
M gas and M sta rs to Mjjm to give a new Ma- A new r& 
is then computed, and the process repeated, until r& 
changes b y < 0.001%. (We note tha t in our previous 
studi es by iHumphrev et alj 120061 and iZappacosta et al.1 
120061 we also computed the virial radius appropriate for 



5 The adiabatic contraction code we used 
was made publicly availa ble by Oleg G nedin at: 



|http://www. astronomy.ohio-state.edu/~ognedin/contra/ 1 
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NGC 5044 




Radius (kpc) Radius (kpc) Radius (kpc) 



NGC 1550 




NGC 2563 




Radius (kpc) Radius (kp C ) Redius (kpc) 



A 262 




Fig. 3. — Results for the emission-weighted projection of the gas density squared (left panels), the emission-weighted projected temperature 
(central panel), and the total gravitating mass (right panel) for NGC 5044, NGC 1550, NGC 2563 and A262. In the temperature and density 
plots, red symbols corresponds to Chandra data, while black symbols corresponds to XMM data. Residuals from the best-fit "parametric 
mass method" models (SjSJl for NFW(+stars) are also shown. In the gravitating mass plot the different mass components are shown: DM 
with the dotted (blue) line, gas mass with the dashed (green) line and stellar mass with the dotted (red) line. Representive "data points" 
are plotted in the gravitating mass profile to show the size of the error bars on the total gravitating mass. 
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Fig. 4.— Same as in Fig|3]for NGC 533, MKW 4, IC 1860 and NGC 5129. The crossed values for the annular bin around 100 kpc for 
IC 1860 indicates that the data point has not been considered in the fit. 
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7. RESULTS 
7.1. Mass-fitting results 

We tested the different mass models listed in § 16.31 
against the data. In the following analysis we obtain 
the best fit by minimizing x 2 ■ Although our best fit 
models are not formally acceptable, the major contribu- 
tions to x stem generally from the inner data points 
(< 10 kpc), where the errors in both temperature and 
density are small. It is not expected that the DM halo 
of each s ystem will be perfectly fitted by an NFW pro- 
file (e.g.. iTasitsiomi et al.ll2004l ). Consequently, we also 
quote the values of the maximum fractional deviation 
dvimax, which gives equal weight to all radial bins, as 
a figure of goodness of fit, in addition to x 2 m Table 
[3] The quantity dvi max is routinely used i n the fits to 
halos formed in numerical simulations; e.g.. I Jin j (2000) 
proposes that dvi max < 0.30 represents a good fit of 
the NFW model. The results for the best fit NFW or 
NFW+stars model are listed in Table [4] at the appropri- 
ate over-density covered by the data. 

Our basic result is that the NFW model, sometimes 
benefiting from an additional component from the stellar 
mass in the central galaxy, is a good overall description of 
the mass profiles. While the formal quality of the fits, as 
noted above, are generally not acceptable, the fractional 
deviations of the fits are typically ~ 10%. The largest 
maximum deviation is observed for NGC 5044 within its 
central radial bin (< 3 kpc), where the Chandra ima- 
ge shows i rregularities presu mably associated with AGN 
feedback (|Buote et al.1 12003). At all other radii the de- 
viations are < 10% for NGC 5044. 

The stellar mass component is not uniformly required. 
When using the NFW+stars model, only 8 objects of the 
16 in the sample show an improvement in the fit. The 
improvement of the fit is judged by considering a reduc- 
tion in x 2 an d a reduction in fractional residuals. This 
provides a quantitative assessment of the improvement 
of the fit even if the final % 2 is still not formally accept- 
able. For example an NFW fit to the MKW4 mass profile 
gives x 2 /dof=58/25 with a dvi max of 0.60 arising from 
the central density bin, while the best fit NFW+stars 
gives x 2 /dof=34/24 with a dvi max of 0.13, because the 
inner data points are better modeled. 

Moreover the amount of improvement is sensitive both 
to the number of data points sampling the inner sa 20 
kpc (where the stellar mass is expected to make a sub- 
stantial contribution to the total mass budget) and to the 
luminosity of the central galaxy. For this purpose it is 
instructive to examine those systems that require stellar 
mass and have both Chandra and XMM data — NGC 
1550, A 262, NGC 533 and MKW 4. By fitting only 
the XMM data, with its coarser binning at small radius, 
we can assess the importance of having high-resolution 
Chandra data for detecting a stellar mass component. 
When fitting only the XMM data for these systems the 
evidence for a stellar mass component is weaker, and the 
inferred amount of stellar mass less, than for the simul- 
taneous Chandra-XMM fits. The amount of stellar mass 
inferred is always larger when the Chandra data are in- 
cluded. In 3 of the 4 cases, the derived concentration 
value does not change within the 1 — 2a errors. The ex- 
ception is A262 for which ca=1-2±0.1 is obtained using 
only XMM and ca=2. 1+0.2 is found for the simultaneous 



Chandra-XMM fits. 

It follows that for systems having only XMM data it is 
necessary to obtain high-quality Chandra observations to 
make a reliable detection of stellar mass. There is clear 
failure to detect stellar mass in 3 objects in our sample 
that are adequately sampled by Chandra observations: 
NGC 5044, NGC 4325 and RGH 80. 

The omission of the stellar compo nent in the mass 
fits leads to biased high conce ntrations (jMamon fc Lokasl 
120051: iHumphrev et al.ll2006f ). but the relevance of the 
bias depends on the number of data points sampling both 
the stellar component (dominant in the inner ~ 20 kpc) 
and DM component. The objects in our sample have 
adequate sampling of the DM component at relatively 
large radii, but the stellar component is well sampled 
(~3 data bins in the inner 20 kpc) only when Chandra 
data are present. As a consequence, the bias is more pro- 
nounced when Chandra data are included. This effect is 
most evident for MKW4. An NFW fit to the XMM data 
for MKW4 gives c A = 5.8 ± 0.3 while an NFW+stars 
fit gives ca = 4.8 ± 0.4. If we use Chandra and XMM 
data, then the fit is driven by the increased number of 
data points within 20 kpc. Fitting an NFW model yields, 
ca = 6.8 + 0.2, which represents a 58% increase over our 
best fit NFW+stars value, c A = 4.3 ± 0.3 (see Table |4j). 
For the remaining objects, fitting only the NFW model, 
when NFW+stars is required, returns a ca biased high 
in the range 38% (A262) to 10% (NGC 533). As ex- 
pected, the bias is generally less for our groups-clusters 
(M > 10 13 M©) than obtained for the e lliptical galaxies- 
group s (M < 10 13 Mq) analyzed by IHumphrev et al.l 

pool) . 

In order to explore the presence and relevance of adi- 
abatic contraction we fitted an NFW*AC+stars model 
to the 8 objects which require a stellar mass compo- 
nent, because only for these objects is the AC model 
potentially relevant. The quality of the fits is not im- 
proved by the introduction of AC (see Table [3]). Be- 
cause the AC model increase the cuspiness of the DM 
profile, we find that the stellar mass (and the derived 
stellar mass-to-light ratios, see Table [5]) were consider- 
ably lower for the NFW*AC+stars model than for the 
NFW+stars model. Because less stellar mass is obtained 
for the AC models, the derived ca values increase by 
10-40% compared to NFW+stars. Two exceptions are 
MKW 4 and RXJ 1159.8+5531, for which AC increases 
ca to 7.1 ± 0.4 and 9.6 ±1.9 respectively. The quality of 
the NFW*AC+stars fits is considerably worse in these 
two cases compared to NFW+stars. The Ma values ob- 
tained for the NFW*AC+stars model are also lower by 
5-20%, with a maximum of 33% for RXJ 1159.8+553. 

Finally, we examined the N04 model. We explored 
N04+stars because in our previous a nalysis of the cluster 
Abell 2589 (jZappacosta et al.ll2006h N04 allowed for an 
increased contribution from stellar central mass compo- 
nents (with values of the Sersic parameter a ~ 0.4). Even 
if we left the Sersic parameter a free, the fit improved 
only in few cases - and only in two, A 2717 (x 2 /dof=7/4) 
and IC 1860 (x 2 /dof=7/4) was the improvement superior 
to 90% according to the F-test. The inferred values of a 
for the sample were quite large and incomp atible with the 
mean value of 0.17 ± 0.03 for CMD halos (jNavarro et al.l 
120041) . If we fixed the value of a at 0.17 the fits did not 
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Same as in Fig|3]for NGC 4325, ESO 5520200, AWM 4 and ESO 3060170. 
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improve. 

7.2. Stellar Mass-to-Light ratios 

Using the stellar mass derived from our fits we cal- 
culated the stellar M/L ratios (M*/L) for the cen- 
tral galaxy. The optical lumin osities have been calcu - 
late d in the Ka-band, following iKochanek et al.l (|2001l ) 
and iLin k Mohrl (|2004l ). using (1) the 20 mag arcsec -2 
isophotal elliptical aperture magnitude, (2) the value 
of the Galactic extinction provided by NED, (3) a k- 
correction of the form k{z) = — 61og(l + z), and (4) 
a correction of 0.2 mag to convert between the total 
and isophotal absolu te magnitudes (see Appendix of 
IKochanek etHl l200lh . We compare this estimate to 
the total extrapolated magnitudes listed in the Extended 
Catalog, finding agreement to better than 10% 6 . For dis- 
tances we adopted the luminosity distance listed in Table 
[TJ Magnitudes have been converted to uni ts of B and 
Ks so lar luminosities using Mb© = 5.48 (|Girardi et al.l 
l2000j) and MksQ — 3.34, which follow s from adopting 
(B- V)Q = 0.64 and (V-Ks)Q = 1.50 dHolmberg et all 
120011 . 

We use the K-band to quote M^/L because Near- 
infrared (NIR) luminosities are much more closely cor- 
relat ed with the total gal axy mass than optical luminosi- 
ties (|Gavazzi et al.|[l996l ). Table [5] shows the best fitting 
results for NFW+stars and NFW*AC+stars for those 
objects requiring a stellar mass component in the mass 
analysis in £17.11 

8. RESULTS FOR INDIVIDUAL OBJECTS 

In the right panels of Fig[3l FigEJ Figj5]and Fig|6]we 
show the total gravitating mass profiles for the objects in 
our sample, with the different components (DM, gas mass 
and stellar mass of the central galaxy) shown in different 
colors and line styles. Details for some individual objects, 
and comparison with previous results in the literature, 
are provided below. 

NGC 1550. Our density an d temperature pr ofiles agree 
with the Chandra analysis of I Sun et~a l. (2003). Their fit 
to the total mass profile within 200 kpc, not surprisingly, 
prefers a Moore profile over an NFW (in particular in 
the inner 10 kpc) because of the stellar mass contribution 
which steepens the profile of the total gravitating matter. 
Our derived DM parameters for an NFW fit are not very 
different from their NFW fit to the total matter: our 
scale radius r s =48 ± 4 kpc and 8 C = 7.76 ± 0.56 x 10 4 
are similar to their best-fitting values of 41.8 kpc and 
1.10 x 10 5 respectively. 

Abell262. The total mass profile has been ana- 
lyzed us ing both Chandra and ROSAT data at large 
radii by IVikhlinin et~aT1 (|2006h . who find a concentra- 
tion, C500 = 3.54 ± 0.30 which is consistent within 
2cr with ours (4.5 ± 0.4, see Table [7]). Other rele- 
vant quantities are in excellent agreement: our values 
of M 2 5oo=3.59±0.14 x 10 13 M Q , / SO s,2500= 0.072 ±0.001 
and r 50 o=624 ± 15 kpc agree well with their values of 
3.40 ± 0.50 x 10 13 M 0! 0.067 ± 0.003 and 650 ± 21 kpc. 

s For a discussion regarding the use of the ellipti- 
cal isophotal magnitude instead of the extrapolated total 
magnitude, because it is less vulnerable to stellar con- 
tamination and surface brightness irregularities, see the 
FAQ sheet for the 2MASS Extended source catalog at 



Using XMM data fPiffaretti etldl (|2005f) obtained the fol- 
lowing parameters from a single NFW model fitted to 
the gravitating matter: r s =85 ± 17 kpc, 0200=8.6 ± 1.0, 
Mfl M , 25 oo=l-97±0.27 x 10 13 M and M ffOS , 25 oo = l-36 ± 
0.20 x 10 12 M . These values do not agree with ours, 
even when we similarly fit only the NFW model (i.e., no 
separate accounting for stars or gas) to the XMM data: 
r s =174± 10 kpc, -Mdm,2500=3- 39 ±0.10 x 10 13 M and 
Af ffQS! 25oo=2.76±0.06x 10 12 M Q . Consi dering the agree- 



ment between our results and those of [ Vikhlinin et al.l 
(120061) . it is unclear why iPiffaretti et al l (|2005h obtain 

different results for this system. 

NGC 533. Using XMM data lPiffaretti et all ([2001 ob- 
tain the following results for the NFW model fitted to the 
gravitating mass: r s =37 ± 3 kpc and C2oo=12.53 ± 0.55. 
Under the same conditions we obtain good agreement 
with their results: r s =43 ± 4 kpc and C2oo=13.0 ± 0.9. 
MKW 4- Our derived temperature pro file shows a de- 



clining behavior like the one obtained by Vikhlini n et al 
(2005), using Chandra data, by iFukazawa et al 



|http://spider.ipac.caltechTe du / staff /j arrett /2mass /XSC/j arrett_XSC primer.html| 



(2004 ), using bo t h Ch andra and XMM data and by 
iFinoguenov et alj (|2007t ) using XMM data. However, our 
profile does not agree with the relativel y flat profile with 
higher temperature values obtained by lO'Sullivan et all 
( 2003]). We believe the origin of the discrepancy proba- 
bly arises from their application of the double subtrac- 
tion method to subtract the background (£ i3.3p . because 
the emission from MKW 4 fills the entire XMM field 
of view. The subtraction of a source component artifi- 
cially hardens the ou termost annuli. A simila r conclusion 
has been reached bv IFinoguenov et"aT1 (|2007[ K Our mass 
model extrapolated to A = 500 gives, C500 = 6.4 ± 0.5, 
r 50 o = 527 ± 8 kpc and A/500 = 4.27 ± 0.18 x 10 13 
M n, which do no t agre e with the parameters found 
bv IVikhlinin et alj (|2006D . c 500 = 2.54 ± 0.15, r 500 = 
634 ± 28 kpc and M 500 = 7.7 ± 1.0 x 10 13 M obtained 
by combining Chandra and ROSAT data at large radii 
(out to 550 kpc). However, when restricting the com- 
parison to the radial range covered by our data, our 
mass (A/2500 = 2.4 ± 0.1 x 10 13 M Q ) agrees with theirs 
(A/2500 = 2.8 ± 0.3 x 10 13 M ). The difference in con- 
centrations stem primarily from a difference between our 
measured scale radius, r s = 81 ± 7 kpc and their value 
of 250 kpc. As we discuss in 19.71 a measurement of 
the scale radius is reliable only if it is well within the 
radia l range of the d ata. Although by using ROSAT 
data ["Vikhlinin et al.l (|2006f ) have surface brightness in- 
formation out to 550 kpc, accurate spectral information 
is available only with Chandra data, which beyond ~ 100 
kpc (outside of the ACIS-S3 chip) are only covering a sec- 
tor of the entire radial annulus and have relatively low 
S/N. Th e similar scale radius of 7 6 kpc for this object de- 
rived by iRines fc Diaferiol (2006) using galaxy redshifts 
and identifying caustics in redshift space supports our re- 
sults j_j!Hioji£htiievalueof r 500 = 620 kpc agrees better 
with ! Vikhlinin et all (|2006l ) (634 ±28 kpc) than with the 
present work (527 ± 8 kpc). Other possible reasons for 
the discrepancy are the different mass modeling proce- 
dure or the different radial range used in the fit, r e strict ed 
to r > 37 kpc in the analysis of Vikhlinin et al. (2006). 

IC 1860. The group exhibits a sharp decline in surface 
brightness in the NW and enhanced emission in the SE. 
This particularly affects the annulus between 94 and 121 



Dark Matter and Gas Fraction in Galaxy Groups. 



15 



TABLE 3 
Quality of the mass fits 



Group x 2 dvi max 

NFW(+stars) NFW*AC+stars NFW(+stars) 



NGC 5044 228/20 0.83 

NGC 1550 66/33 74/33 0.22 

NGC 2563 18/7 0.23 

A 262 116/37 134/37 0.19 

NGC 533 81/20 85/20 0.30 

MKW 4 34/24 63/24 0.13 

IC 1860 14/5 25/5 0.11 

NGC 5129 3/2 3/2 0.15 

NGC 4325 12/9 0.29 

ESO 5520200 18/8 0.08 

AWM 4 13/9 14/9 0.11 

ESO 3060170 25/9 0.12 

RGH 80 71/9 0.42 

MS 0116.3-0115 6/5 0.16 

A 2717 32/9 0.22 

RXJ 1159.8+5531 5/7 13/7 0.17 



Note. — Various indicators of quality of fit for the mass models discussed 
in the text. 

X 2 refers to the x 2 /dof of the fits to the density and temperature profiles used 
to infer the parameters of the mass model in the gas potential approach. 
dvi max refers to the maximum fractional deviation on the fits to the density 
and temperature profiles used to infer the parameters of the mass model in 
the gas potential approach. 

Values of x 2 for the NFW*AC+stars model are reported only for the objects 
showing an excess over the NFW fit due to stellar mass. 



TABLE 4 

Results for the NFW virial quantities at selected overdensity 



Group 


A 


r a 


CA 








f gas, A 


Mdm,a 








(kpc) 




(kpc) 


(10 13 M ) 


(10 12 M Q ) 




(10 13 M Q ) 


(10 10 M ) 



NGC 5044 


1250 


77 + 2 


3 


.8 + 0. 


1 


295 ± 2 


1.85 + 0.04 


1 


.21 + 0.02 





.065 ± 0.001 


1 


72 ± 0. 


04 




NGC 1550 


2500 


48 + 4 


4. 


.5 + 


.3 


215 ± 2 


1.42 + 0.03 


1 


.02 + 0.02 





072 ± 0.001 


1 


31 + 


03 


11.2 + 4.1 


NGC 2563 


2500 


76 + 22 


2. 


.4 + 1 


.0 


185 ± 5 


0.92 + 0.08 





.31 + 0.03 





.034 ± 0.001 


0. 


89 + 0. 


08 




Abell 262 


2500 


141 ± 16 


2. 


.1 + 


.2 


292 + 4 


3.59 + 0.14 


2 


.60 ± 0.08 





.072 ± 0.001 


3 


.31 ± 0. 


.13 


22.1 + 4.5 


NGC 533 


1250 


43 + 4 


6. 


.1 + 


.5 


262 ± 2 


1.30 + 0.04 





.87 + 0.02 





.067 ± 0.001 


1 


.19 + 0. 


.04 


22.4 + 2.2 


MKW4 


1250 


81 + 7 


4. 


.3 + 


.3 


353 ± 4 


3.21 + 0.10 


2 


.84 ± 0.06 





.088 ± 0.002 


2 


.87 + 


.10 


61.8 + 7.2 


IC 1860 


1250 


101 ± 12 


3 


.2 + 


.3 


319 ± 6 


2.36 + 0.13 


1 


.56 + 0.05 





.066 ± 0.002 


2 


.18 + 


.12 


26.4 + 6.3 


NGC 5129 


1250 


43 ± 10 


5. 


.2 + 


.9 


226 ± 7 


0.84 + 0.07 





.58 + 0.06 





.069 ± 0.003 





.78 + 


.07 


9 S+6.7 


NGC 4325 


2500 


75 ± 18 


2. 


.8 + 


.4 


208 ± 8 


1.32 + 0.16 





.66 ± 0.03 





.050 ± 0.004 


1 


.26 + 


.16 




ESO5526020 


1250 


171 ± 27 


2. 


.5 + 


.3 


422 ± 13 


5.51 + 0.51 


3 


.35 + 0.18 





.061 ± 0.002 


5 


.17 + 


.50 




AWM 4 


1250 


154 ± 17 


3 


.0 + 


.3 


465 ± 13 


7.38 + 0.61 


4. 


.79 + 0.29 





.065 ± 0.003 


6 


.88 + 


.59 




ESO3060170 


2500 


162 ± 54 


2. 


.1 + 


.3 


343 ± 18 


5.97+1.14 


3 


.45 + 0.17 





.058 ± 0.005 


5 


.62 ± 1 


.12 




RGH 80 


500 


78 + 8 


5. 


.1 + 


.5 


397 + 5 


1.85 + 0.07 


2 


.85 + 0.11 





154 + 0.003 


1 


.56 + 


.06 




MS 0116.3-0115 


1250 


202 ±115 


2. 


.0 + 


.8 


405 ± 42 


4.92 ± 1.64 


1 


.97 + 0.19 





.040 ± 0.009 


4. 


.73 ± 1 


63 




Abell 2717 


500 


233 ± 18 


3 


.0 + 


.2 


710 ± 11 


10.68 + 0.51 


11.36 + 0.29 





106 ± 0.003 


9 


.55 + 


.49 




RXJ 1159.8+5531 


500 


104 ± 77 


5. 


.6+1 


.5 


584 ± 73 


6.13 + 3.30 


5.10 + 0.41 





.083 ± 0.019 


5 


.57 + 3 


.32 


56.9 + 10.5 



Note. — Results for the mass profile fits. A refers to the overdensity chosen for the object, as the closest to the outer radius of the data. r s is 
the scale radius of the NFW profile. 



kpc, which has been excluded from the fit. We studied 
the effects of this asymmetry by dividing the annuli into 
two sectors. We defined the SE sector as 15-195 degrees 
measured from the N. The corresponding NW sector is 
then 195-15 degrees. We find that the gas density profile 
is steeper in the NW direction, but the lower tempera- 
ture in the 91-125 kpc annulus is caused by the cooler, 
higher density emission in the SE. The radial tempera- 
ture profile is quite smooth over the NW sector. The Ca 
and Ma values obtained for each sector when excluding 



the 94-121 kpc annulus are consistent within their ~ 2a 
errors. Including this annulus has negligible impact on 
the results for the NW sector. These low level distur- 
bances did not indicate a significant violation of hydro- 
static equilibrium, as further suggested by the agreement 
of the derived c v i r and M v - lv with the values expected from 
ACDM simulation. 

NGC 4325. We measured an NFW scale radius r s — 
75 ± 18 kpc and M 200 = 3.01 ± 0.65 x 10 13 . The 
results agree with the uncertain values obtained by 
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TABLE 5 
Stellar mass-to-light ratios 



Galaxy L K /L B Fitted M*/L K (Mq/Lq) 

NFW+H90 NFW*AC+H90 

NGC1550 4.8 0.53 ± 0.20 0.24 ±0.01 

NGC708 10.7 0.54 ±0.11 0.14 ±0.10 

NGC533 4.9 0.36 ± 0.03 0.11 ±0.01 

NGC4073 4.6 0.86 ±0.10 0.33 ± 0.04 

IC1860 7.2 0.60 ±0.14 0.26 ± 0.01 

NGC5129 4.1 0-Oeto.ll 0.05 ± 0.01 

NGC 6051 7.6 0.301q|o 0.18 ±0.01 

2MASSX J11595215 4.4 0.55 ± 0.10 0.40 ± 0.05 



Note. — K-band stellar mass-to-light ratios for the central galaxy 
measured from our fits to the data using the NFW+stars and the 
NFW*AC+stars models. 



iRines fc Diaferiol (|2006l ); i.e., M 200 = 1.5 ± 1.3 x 10 13 
and r s — 82 kpc. 

AWM4- This object has a remarkable temperature 
profile. Unlike the other groups in our sample, the core 
is isot hermal as found previously by 10 'Sullivan et al.l 
(2005). Beyond a radius of 200 kpc' we measure a de- 
clining temperature profile with the XMM data. We 
find that the source emission fills the entire field o f 
view, contrary to the analysis in lO'Sullivan et all ((2005) , 
which reported a "soft excess" described by a 0.6 keV 
bremsstrahlung component, probably the misinterpreted 
source. It is difficult to classify AWM 4 as a merging 
system, given its r elaxed appearance both in the X-rays 
and in the optical (jKoranvi fe Gellerll2002l ) . Instead, the 
flat temperature profile likely reflects the influence of the 
powerful AGN, with radio lobes extending out from the 
central galaxy NGC 60 51 along the minor a xis of the 
galaxy to 100 kpc (e.g.. iNeumann et al.lll994D . 

ESO 3060170. Our temperature profile is best fit- 
ted by a declining profile at large radii. However, be- 
cause of the relatively large error bars, our profile is 
also co nsistent with the flat profile obtained bv lSun et al.l 
(120041) between 10 and 400 kpc with XMM and Chandra. 
iSun et alJ (|2004[ ) obtain c 200 - 8.7 and M 2Q0 ~ 1.8 x 10 13 
M Q which may be compared to our extrapolated values, 
caoo = 6.7 ± 0.8 and M 200 = 1.54 ± 0.59 x 10 14 M Q . 

RGH80. The Chandra image clearly reveals the peak 
of the X-ray emission coincident with MCG ±06-29-077 
and a bright tail pointing NW - with MCG ±06-29-078 at 
the S edge of this feature. This geometry was only hinted 
at by the XMM image (see Fig. 10 of iMahdavi et all 
2005). This asymmetry is an indication that the core is 
not fully relaxed, as already suggested by the absence 
of a single central galaxy. Despite this fact, hydrostatic 
equilibrium beyond the inner core seems a good approx- 
imation given the values of c v i r and M v j r measured for 
this object. 

Abell2717. The temperature profile we have derived 
from XMM data declines at large radii like all of the 
groups in our sample and is inc o nsiste nt with the flat pro- 
file found bv iPratt fc Arnaudl §005). The origin of the 
difference is likely our improved treatment of background 
subtraction. Nevertheless, our inferred c 2 oo = 4.6 ± 0.2, 
M200 = 1-59 ± 0.06 x 10 14 Mq and r 200 = 1082 ± 21 
are in good agreement wit h those determined both b y 
IPratt fc Arnaudl (|2005l ) and iPointecouteau et all (|2005h . 
C200 = 4.21 ± 0.25, M 200 = 1.57 ± 0.19 x 10 14 M and 



^200 = 1096 ± 44. The reason for the agreement, despite 
the difference in the temperature pr o files, is likely the 
same put forward bv IVikhlinin et all (|2006h : the NFW 
fit implies a declining temperature profile at large radii. 

RXJ 1159+5531. Our i nferred c 5 qq = 5 6 ± 1.5 is higher 
than the one reported in IVikhlinin et alJ (|2006l ). C500 = 
1.7 ± 0.3, using the same Chandra data, though within 
2.5cr given our large error bars. Our derived Af 2 5oo = 
3.3 ± 0.9 x 10 13 Mq and gas fraction f gas ,2500 = 0.049 ± 
0.004 are on the contrary in good agreement with their 
determination of 3.0±0.3 x 10 13 M Q and 0.045±0.002. As 
for MKW 4 the key difference is in the measured scale 
radius: our value of 104 ± 77, though with large error 
bars, is inconsistent with their quite high value of 412 
kpc, which is again at the boundary of the radial range 
covered by the data, which are of not excellent quality 
outside the S3 chip (~ 370 kpc). 

9. SYSTEMATIC ERRORS 

In this section we address the sensitivity of our analy- 
sis to various systematic uncertainties and data-analysis 
choices which may impact upon our results. An estimate 
of the uncertainty due to these effects for each object 
is given in Table El The statistical error of the default 
model (AStatistical) is also listed on the table. In the 
case of the different approach of using a deprojection tec- 
nique ( ADeproj), we also quote the corresponding ma- 
gnitude of statistical error together with the associated 
best-fitting parameter shift in the table. We illustrate 
the effect of systematic errors on the best fit ca, Ma 
parameters, and the stellar mass-to-light ratio M±/Lk- 

9.1. Background modeling and subtraction 

One of the most important potential sources of sys- 
tematic uncertainties in measuring the mass profiles of 
groups is the background subtraction technique, in par- 
ticular in the low surface brightness regime at large radii. 
Our modeling technique is particularly effective in the 
low temperature regime of groups, and we take as an 
extreme measure to change the overall normalization of 
the background model by ±5%. Such an error in the 
estimated background is unlikely, but the exercise is in- 
dicative of our sensitivity to the background. 

9.2. Spectral-fitting choices 

Among the variety of choices made in spectral-fitting, 
we explore the ones more likely to affect to some degree 
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TABLE 6 
Systematic Error Budget 



Group 


Best Fit 


AStatistical 


ABackground 


ASpectral 


AMethod 


ADcproj 


Ar e 






NGC 5044 


3.8 


±0.1 


-0.5 


±0.2 


±0.6 


-0.3 (±0.2) 




NGC 1550 


4.5 


±0.3 


±0.4 


±0.2 


-0.1 


-0.8 (±0.3) 


±0.5 


NGC 2563 


2.4 


±1.0 


±2.6 


±2.3 


-0.1 


+4.5 (±1.4) 




Abell 262 


2.1 


±0.2 


±0.2 


+0.8 
-0.6 


-0.4 


-0.4 (±0.2) 


-0.2 


NGC 533 


6.1 


±0.5 


-1.7 


-2.0 


±1.1 


-1.5 (±0.4) 


±0.9 


MKW 4 


4.3 


±0.3 


-0.1 


+0.3 
-0.7 


-0.3 


+0.8 (±0.7) 


-0.3 


IC 1860 


3.2 


±0.3 


+0.1 


+0.9 
-0.4 


-1.3 






NGC 5129 


5.2 


±0.9 


±0.6 


-0.4 


-0.3 


-0.2 (±2.2) 




NGC 4325 


2.8 


±0.4 


±0.7 


±0.9 


±0.3 


-0.7 (±0.3) 




ESO5520200 


2.5 


±0.3 


-0.2 


-0.3 


±0.1 


+0.2 (±0.4) 




AWM 4 


3.0 


±0.3 


+0.1 


-0.2 


-0.1 


-0.9 (±0.3) 




ESO3060170 


2.1 


±0.3 


-0.4 


+0.8 
-0.6 


-0.3 


-0.1 (±0.3) 




RGH80 


5.1 


±0.5 


±2.1 


±4.5 


-2.6 


+2.9 (±1.2) 




MS 0116.3-0115 


2.0 


±0.8 


+0.7 


+ 1.5 
-0.5 


±1.0 


+2.3 (±1.9) 




Abell 2717 


3.0 


±0.2 


+0.1 


-0.2 


-0.1 


+0.6 (±0.3) 




RXJ 1159.8+5531 


5.6 


±1.5 


-0.9 


±0.7 


-1.2 


+2.6 (±1.7) 




M A /10 13 M Q 


NGC 5044 


1.85 


±0.04 


±0.28 


-0.10 


-0.41 


±0.34 (±0.09) 




NGC 1550 


1.42 


±0.03 


-0.04 


-0.03 


±0.02 


±0.26 (±0.09) 


±0.01 


NGC 2563 


0.92 


±0.08 


-0.06 


-0.17 


±0.01 


-0.24 (±0.13) 




Abell 262 


3.59 


±0.14 


-0.19 


+0.24 
-0.62 


±0.34 


±1.00 (±0.31) 


±0.10 


NGC 533 


1.30 


±0.04 


+0.15 
-0.01 


±0.16 


-0.04 


-0.01 (±0.07) 


-0.05 


MKW 4 


3.21 


±0.10 


-0.10 


+0.12 
-0.07 


±0.09 


-0.86 (±0.18) 


±0.03 


IC 1860 


2.36 


±0.13 


-0.08 


+0.12 
-0.20 


±0.65 






NGC 5129 


0.84 


±0.07 


+0.08 
-0.03 


-0.02 




-0.13 (±0.15) 




NGC 4325 


1.32 


±0.16 


-0.15 


-0.20 


-0.10 


±0.53 (±0.45) 




ESO5520200 


5.51 


±0.51 


±0.35 


+0.70 
-0.13 


-0.40 


±0.49 (±0.71) 




AWM 4 


7.38 


±0.61 


-0.27 


-0.70 


±0.16 


±2.01 (±0.87) 




ESO3060170 


5.97 


±1.14 


±1.30 


+2.07 
-0.74 


±0.73 


±0.68 (±1.37) 




RGH80 


1.85 


±0.07 


+0.26 
-0.14 


+0.05 
-0.40 


±0.48 


-0.07 (±0.19) 




MS 0116.3-0115 


4.92 


±1.64 


±0.46 


+ 0.68 
-1.42 


-1.12 


-0.40 (±3.76) 




Abell 2717 


10.68 


±0.51 


-0.03 


±1.02 


±0.49 


-0.76 (±0.86) 




RXJ 1159.8+5531 


6.13 


±3.30 


±0.97 


-0.29 


±0.51 


-1.87 (±0.72) 





M,/L K (M Q /L Q ) (NFW±stars) 



NGC 1550 


0. 


.53 


±0.20 


±0.05 


±0.12 


±0.03 


±0.17 (±0.15) 


±0.36 


Abell 262 





.54 


±0.11 


-0.13 


+0.23 
-0.37 


+0.12 


-0.45 (±0.04) 


±0.79 


NGC 533 





.36 


±0.03 


±0.20 


±0.26 


-0.27 


±0.24 (±0.04) 


-0.07 


MKW 4 





.86 


±0.10 


+0.51 
-0.66 


+0.59 
-0.54 


+0.11 


±0.15 (±0.09) 


±0.60 


IC 1860 


0. 


.60 


±0.14 


±0.01 


+0.02 
-0.23 


±0.53 






NGC 5129 





.06 


+0.13 
-0.06 


+0.13 
-0.06 


-0.06 


±0.06 


±0.43 (±0.19) 




AWM 4 





.30 


+0.33 
-0.30 


±0.08 


±0.12 


±0.14 


±0.82 (±0.38) 




RXJ 1159.8+5531 





.55 


±0.10 


±0.05 


+0.10 
-0.13 


+0.11 


-0.13 (±0.13) 





Note. — The estimated error budget for each of the groups. Excepting the statistical error, these values estimate a 
likely upper bound on the sensitivity of the (best fit) value of each parameter to various data-analysis choices and should 
not be added in quadrature with the statistical error. The "Best" column indicates the best-fit value and "AStatistical" 
the lcr statistical error for Ma and ca from Tabic |4| and for M„/L K from Tabled "ABack ground" gives the results when 
the X-ray background level is set to ±5% of nominal. "ASpectral" gives the results when changing spectral analysis 
choices, AMethod when adopting a different approach (using Eq[2]or Eq[3j to mass modeling, ADcproj when using 
projected (2D) or dcprojcctcd (3D) spectral results (with the associated statistical error), and Ar e when changing the 
effective radius of the stellar profile. 
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Same as in Fig|3]for RGH 80, MS 0116.3-0115, A 2717 and RXJ 1159.8+5531. 
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the inferred gas density and temperature in each radial 
bin. 

The plasma code. Different plasma codes choose from 
a large, overlapping, but incomplete set of atomic data, 
leading to differences in the inferred abundances and, 
therefore, density and, to a lesser extent, temperature. 
We experimented with replacing the APEC model with 
the MEKAL plasma model. 

Bandwidth. To estimate the impact of the bandwidth 
on our fits, we experimented with fitting the data with 
different lower limits for the energy band. In addition to 
our preferred choice of 0.5 keV, we use 0.4 keV and 0.7 
keV. 

Hydrogen column- density. We take into ac- 
count possible dev i ations for N-r from the value of 
iDickev &: Lockmanl (|1990f ) allowing the parameter to 
vary by ±25%. 

9.3. Deprojection method 

We analyzed the possible systematics involved with 
the projection of 3D mo dels using instead the "onion- 
peeling" technique (e.g., iFabian et alJll98ll iKriss et all 
[1981 iBnotejl2000bft . Only for the object IC 1860 we 
did not perform this exercise because of the exclusion of 
an inner bin (see $8|). The results were consistent with 
the ones obtained by the 2D analysis (see Table [6|) but 
with larger error bars given the quality of the current 
data. This is the main reason for having adopted the 2D 
analysis as our default. In Fig|7] we plot as a function 
of the fraction of the virial radius the quantities Ap/p^jj 
(where Ap = p2D — Pad with p^o the value of the best fit 
model of the 2D analysis) and the corresponding quan- 
tity AT/T3D. The plotted errors are the fractional errors 
on the derived 3D quantities. There is more scatter in 
the density as a consequence of larger uncertainties in 
the derived 3D iron abundances, while the temperatures 
determined with the two methods are generally consis- 
tent. This fact reinforces the notion that, for the range of 
temperatures spanned by the objects considered in this 
paper, the spectroscopic temperatures are not biased sig- 
nificantly (see discussion in £19.41 ^10.5|) . 

9.4. Response weighting 

Since the effective area of the detector response of 
ACIS on Chandra and EPIC on XMM are a decreas- 
ing function of energy, fitting a IT model to a spec- 
trum having a range of temperature components above 
~ 3 keV will tend to yield a temperature that is biased 
low (jMazzotta et al.l l2004t IVikhlininl 120061) . This effect 
is negligible for average temp eratures around 1 keV (see 
Appendix [Bland Buotc 2000c) as it is the case for the ob- 
jects in our sample. As a further systematic check we ap- 
plied a straightforward averaging of the plasma emissiv- 
ity over the detector response in our method as explained 
in Appendix [Bl The results obtained using the response 
weighting are very consistent with the ones obtained by 
our default projection analysis. For example, for NGC 
1550 we obtain c A = 4.5 ± 0.3, M A = 1.45 ± 0.03 x 10 13 
M Q and M*,a = 11.1 ± 4.0 x 10 11 M Q ; for MKW 4 we 
obtain c A = 5.1 ± 0.4, M A = 2.91 ± 0.10 x 10 13 M Q 
and M*, A = 60.7 ± 7.3 x 10 11 M Q ; for NGC 533 we ob- 
tain c A '= 5.3 ± 0.4, M A = 1.41 ± 0.06 x 10 13 M© and 
M* )A = 31.6 ±2.1 x 10 11 Mq. 



9.5. Mass derivation method 

For each system we tried all the three methods de- 
scribed in ££6j By using all the approaches we have an 
estimate of the robustness of the inferred mass and virial 
quantities. We also include in this estimate the fact that 
different temperature and density profiles may be able 
to fit the same data adequately but give rise to different 
global halo parameters. To test this, we cycled through 
each of our adopted gas density and temperature profiles. 

9.6. X-ray asymmetries and disturbances 

There are systems displaying low level asymmetries (IC 
1860), substructure (RGH 80), and AGN cavities (A 262) 
or possible AGN-induced disturbances (NGC 5044). For 
the objects which have a mild degree of disturbance in 
the core we found that the results obtained excluding 
the disturbed regions agreed with those obtained over 
the entire radial range within the 1-2<t errors. For A262 
we excluded the inner 20 kpc to avoid (1) the cavities 
which affect the central 10 kpc, and (2) the stellar mass 
component of the central galaxy. In this case fitting an 
NFW profile gives, c A = 2.4 ± 0.3 and M A = 3.44 ± 
0.15 x 10 13 M Q . For NGC 5044 we obtain c A = 3.9 ± 
0.1 and M A = 1.83 ± 0.04 x 10 13 M after excluding 
the central 5 kpc where there is evidence of a disturbed 
morphology. We exclude the inner 30 kpc of RGH 80 
and find c A = 6.5 ± 1.0 and M A = 1.78 ± 0.08 x 10 13 
Mq. Finally, for IC 1860 we perform a sector analysis, 
extracting spectra and re-deriving our mass profiles from 
suitably oriented semi-annuli, as detailed in <JS] We found 
consistent results within their ~ 2a errors. Therefore, we 
infer no systematic error associated with including the 
central, mildly disturbed, regions in these systems. 

9.7. Radial Range & Extrapolation 

It is customary to extrapolate mass profiles out to 
the virial radius defined within an over-density A ~ 
100 — 500. This facilitates a consistent comparison to 
theoretical studies which usually quote results in this ra- 
dial range, corresponding to the entire virialized portion 
of the halo. X-ray studies of global scaling relations be- 
tween mass, temperature, and luminosity also prefer to 
use such large virial radii to seek the tightest relations 
between these global quantities. 

However, extrapolating the mass profiles can lead to 
systematic errors in Cy j r , M v j r , and the gas mass/fraction. 
iVikhlinin et "all ((2006) argue that biased extrapolation of 
the gas density profiles is the main reason for the under- 
estimate of gravitational masses and low normalizations 
of the M — T relations found wit h earlier X-ray tele- 
scopes (e.g.. iNevalainen et al.ll2000l ). using a (3 model fit 
for the gas density and a polytropic appr oximation for 
the temperature profile. iRasia et "ail (|2006l ) suggest that 
the same systematic error affects c V i r , in the sense that 
a restricted radial range tends to return a higher c v ; r , in 
the context of the NFW profile, than the value derived 
using data extending out to the virial radius. 

Our procedure for mitigating extrapolation bias as- 
suming the halo follows an NFW profile is as follows. We 
obtain the mass profile within an appropriate r A corre- 
sponding to the outer radius of the X-ray data for each 
group. The values of c A and Ma are extrapolated to 
A sa 101 assuming the NFW profile applies, using the 
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various systems are the same in both plots and they are listed on the plot in the right panel. 



convenient approximation of IHu Sz Kravtsovl (|2003l ) . We 
emphasize that we do not need to extrapolate the models 
for the gas density and temperature to obtain the extrap- 
olated mass parameters in this manner. The most impor- 
tant requirement for self-consistent extrapolation is that 
the NFW scale radius be accurately measured using the 
available X-ray data at smaller radius. Since our princi- 
pal approach for measuring the mass profile (parametric 
mass method, see $3]) guarantees a physical solution of 
the equation of hydrostatic equili brium for an NFW DM 
hal o, unlike the metho ds used by IVikhlinin et al.l (|2006l ) 
and lRasia et al.1 (|2006f ). and our temperature profiles are 
modeled with a more sophisticated approach than the 
polytropic-/? model estimate, we expect more reliable 
measurements of r s . 

The crucial factor for reliable measurement of r s is 
that the true value of r s lies well within the outer ra- 
dius of the X-ray data. We illustrate this effect using 
those objects for which we obtained measurements out 
to A = 500 (A2717, RGH 80 and RXJ 1159.8+5531). 
If we exclude the outer two data points of A2717 then 
the new outer data point corresponds to 320 kpc and 
A = 2276. Fitting the X-ray data over this mass range 
gives a best-fitting value, r s = 338 kpc, uncomfortably 
outside the new radial range of the data and larger than 
inferred using all of the data (r s = 233 kpc). Extrapo- 
lating this profile to A = 500 yields a larger mass and a 
smaller concentration than obtained for all of the data. 
Analogous results are obtained when performing this ex- 
ercise for RGH 80. For RXJ 1159.8+5531 we exclude 
the outer data point so that the new outer radius is 
289 kpc corresponding to A = 2318. In contrast to 
A2717 and RGH 80, when fitting over this smaller ra- 
dial range, we obtain a scale radius 131 ± 76 kpc, still 
well within the outer radius. The derived mass parame- 
ters are, c 250 o = 2.2 + 0.6 and M 250 o = 3.69 + 0.86 x 10 13 
Mg. Extrapolating these values to A = 500 we obtain 
c 500 = 4.7+ 1.1, M 50 a = 7.10 + 3.25 x 10 13 M Q , in excel- 
lent agreement with the results obtained over the whole 
data range presented in Table [U 

This exercise suggests that measurements of ca and 
Ma should be reliable provided the NFW scale radius 
lies well within the outermost radius covered by the data, 
as is the case for all the objects in our sample. Agree- 



ment with the optical determination o f the scale radius 
for th e two objects in common with iRines &: Diaferid 
( 2006) adds further strength to the results (see $5]). 

Unfortunately, extrapolation of the gas mass and gas 
fraction is less reliable. If we extrapolate our models out 
to a virial radius corresponding to A ~ 101 we obtain 
gas fractions consistent with the cosmic value in 12 of 
16 cases. In 4 cases the extrapolated gas fractions ex- 
ceed the cosmic value derived by WMAP, suggesting a 
problem with the extrapolation. All these systems pos- 
sess a flat slope of the gas density profile ( (3 < 0.5) at 
the edge of the data range. This type of behaviour has 
been noted previously by simulations and simple ana- 
lytic models which pointed out how the [3 model overesti- 
mates gas mass (and underestimates gravitational masses 
based on (3 model fits), because it returns a biased low [3 
due to t he restricted range of radii where the fit is per- 
formed dNavarro et al.l Il995l; [Bartelmann k, Steinmetzl 



IBM [Borgani et alJ 12004 iKomatsu fc SeliakH200ifr Tn^ 
deed, IVikhlinin et al.l (|2006f) finds evidence for a steep- 
ening of the gas density slope with radius in clusters. 



9.8. The stellar mass profile of the central galaxy 

To account for the stellar component we adopted a 
De Vaucoleurs model with effective radius being fixed 
to that determined by 2MASS. The derived stellar mass 
is most sensitive to the effective radius. The difference 
in effective radii measured in different optical bands, 
as evident in Table is mainly due to the use of dif- 
ferent fitting ranges/sensitivity (e.g.. iFisher et aLll!995l) 
and to a radial color gradient, reflecting gradients in the 
metallicity or age of the stellar population (e.g., Pahre 
1999). Though the true stellar mass is more reliably 
determined from K-band data, we investigated the sen- 
sitivity of our parameters, in particular the value for 
M*/L, to the choice of r e , by replacing the K-band r e 
for each galaxy with the larger B-band value, listed in 
the Third Reference Catal og of Bright Galaxies (RC3: 
Ide Vaucouleurs et al.lll99ll ). The stellar mass, and con- 
sequently M*/L, increases systematically, with the only 
exception being NGC 533, when using the larger B-band 
effective radius. But, importantly, the concentration and 
mass are not affected, showing that the main conclusions 
of our paper regarding these two quantities are not sen- 
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sitively dependent on the adopted stellar template. 

We studied the possible contribution to the stellar mass 
of non central galaxies within r s using known member 
galaxies with 2MASS photometry as listed in NED. For 
most of the objects they contribute no more than 10% of 
the total light. This i s cons istent with the more general 
result of iLin fe Mohrl (|2004f ) who showed that the ratio 
of BCG-to-total galaxy light decreases with increasing 
cluster mass indicating that 30-50% of the total light 
in galaxies is in the BCG for group-scale systems (10 13 
M©< M200 < 10 14 Mq). Notable exceptions are infact 
the most massive objects in our sample: A262 has a 77% 
and AWM 4 a 39% additional contribution from non cen- 
tral galaxies within r s . 

10. DISCUSSION 

10.1. Cvir--Mvi r relation 

In Fig. [5]we plot the c vil -M vil relation fits to the XMM 
and Chandra data (A « 101, see Table [7]). To obtain 
an empirical description of the relation we fitted a sim- 
ple power-law mod el following the approach described in 
IBuote et al.l (|2006l) . That is, we fitted the data with a 
linear relation of the form lo g (1 + z)c v j T = a log[M Y j r 1 +b 
using the BCES estimator of lAkritas fc Bershadvl ( 19961 ) 
with bootstrap resampling. We obtain a = —0.226 ± 
0.076, implying that the concentration decreases with 
increasing mass at the 3<r level. The previous studies 
of clusters (> 10 14 M Q ) with Chandra and XMM found 
a sw - very consistent with a constant c v j r -My lr re - 
lation ([Pointecouteau et al.ll2005l : IVikhlinin et al.ll2006[ ). 
Therefore, it is the lower mass range, 10 13 — 10 14 M Q , ap- 
propriate for groups that provides crucial evidence that 
c vir decrease s with increasing My 1T as expected in CDM 
models (e.g.. iBullock et al.ll200lT ). 

The best-fitting power-law model is plotted in Fig. [5] 
along with the theoretica l prediction of the A CDM ob- 
tained using the model of IBullock et al.l (|2001h with pa- 
rameters (F = 0.001, K = 3.12) intended to represent 
halos up to masses ~ 10 14 Mq. Also shown is the pre- 
dicted la intrinsic scatter for the ACDM model. For 
M>4x 10 13 M Q the ACDM model is a good represen- 
tation of the X-ray data. For lower masses, the observed 
c v i r -M V i r appear to exceed the prediction. Allowing for 
a w 10% increase in the concentrations predicted by the 
ACDM model for the most relaxed, early forming ha- 
los (|Jindl2000trWechsler et aHl2002t IMaccio' et al J 1 20061) 
helps to bring the model into better agreement with the 
observations. 

We infer an intrinsic scatter, 0.03±0.02, in logm (l+z)c 
for t he empirical power law relation (see IBuote et al.l 
|2006| ) which is considerably l ess than the value of ~ 0.14 
obtained for ACDM halo s (jJingj 120001 : IWechsler et all 
120021 IMaccio' et aT] 120061) . For the most relaxed, early 
forming halos ACD M simulations typical ly find a smaller 
scatt er ~ 0.10 (jWechsler et al.1 120021 : IMaccio' et aT 
120061) . though the most relaxed halos studied by Uine 
(2000) have a scatter of 0.07 (after converting between 
In and log 10 with a factor of 2.3). The small scatter we 
have measured about the power-law relation only agrees 
with CDM simulations if these halos are the most re- 
laxed, early forming systems. 

The data-model comparison has been made for a 
ACDM model with p arameters obtained from the first 
year WMAP results (|Spergel et alj I2003D . in particu- 



lar with r2 m =0.3 and as = 0.9. With the lower val- 
ues favored by the three year data release of WMAP 
([Spergel et all [20061 ) the pr edicted concentratio ns are 
lower (e.g., see discussion in IMaccio' et al.1 [20 061. The 
implications are discussed in IBuote et al.l ( 2006h . 



10.2. The detection of central stellar mass 

A good fit of the NFW profile to the total gravitating 
matter of relaxed, T > 3 keV, massive clusters, with- 
out any significant deviation arising from the central 
stellar mass, appears to be a common feature of X-ray 
studies (iPointecouteau et "alll2005t IVikhlinin et al.l 



IZappacosta et al.ll2006( ) . On the contrary, relaxed bright 
elliptical galaxie s always require a domi nant contribution 
of stellar mass ([Humphrey et all 12006( 1 . The intermedi- 
ate mass range explored here shows a mixed behavior: 
some low-mass, group-scale objects and three poor clus- 
ters (AWM 4, RXJ 1159 and Abell 262) do show evidence 
of stellar mass , while there are examples of objects whose 
gravitating mass profile is described by just NFW. (Note 
that RXJ 1159, A 262 and MKW 4 were also shown to 
have an excess core mas s profile above that in dicated by 
NFW in the analysis of IVikhlinin et al.ll2006l ) 

An important issue emerging in the analysis is how well 
the two key components, the stellar component associ- 
ated with the central galaxy and the DM, are sampled 
by the X-ray data. It is expected that the stellar com- 
ponent is most relevant within the inner 10-20 kpc while 
the DM should dominate the mass budget elsewhere. To 
reveal and measure adequately the stellar mass, enough 
density and temperature data points are required in the 
inner ~ 20 kpc, depending as well on the amount of stel- 
lar mass present (implied by Lk)- The omission of the 
stellar component in the mass modeli ng, proposed as a 
possi ble source of abnormally high c ([Mamon fc Lokasl 

2005) , is certainly a factor for relatively nearby objects 
with data densely sampling the inner dominated stellar 
core but not extending to larg e radii, as for the object s 
analyzed with Chandra data in iHumphrev et al.l ((2006) . 
The effect is less pronounced in the objects analyzed in 
this sample, where the data extend to large enough radii, 
but with comparatively less density of data points in the 
inner 20 kpc, in particular for objects with only XMM 
observations. The presence of data at large radii pre- 
vent to obtain large values of c v i r (> 30) when fitting the 
wrong NFW model to objects which require stellar mass. 

For the objects which require stellar mass and have 
2-3 data bins in the inner 20 kpc (NGC 1550, A262, 
NGC 533, MKW 4) the derived stellar M/L ratios are 
consistent with the range of values fo und in our anal- 
ysis o f a sample of elliptical galaxies ([Humphrey et al.l 

2006) : the (unweighted) mean M/L ratio of these four 
objects is 0.57 ± 0.21 which is consistent within la with 
0. 76 ± 0.24, the mean ste llar M/L ratio of the objects 
in IHumphrev et al.l (|2006[ ). though a 25% difference is 
present. This result reinforces the relevance of the sam- 
pling of the inn e r reg ion: on average the objects in 
IHumphrev et al.l (|2006l ) have ~ 7 data bins in the in- 
ner 20 kpc allowing a more accurate measurement of the 
stellar mass. The measures for objects with low reso- 
lution XMM observations are likely biased low, as we 
determined for objects with both XMM and Chandra ob- 
servations. 

For a single burst stellar population (SSP) with age 
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ranging from 9-13 Gyr and metallicity ranging from 0.5- 
2 solar, the K-band stellar mass-to-light ratio is ex- 
pected to take values in the range 0.86-1.16 for a Kroupa 
(|Kroupalj200ll) IMF and 1.28-1.49 for a Salpeter IMF 
( Salpeteii Tl955). IMF (We have used linearly interpo- 
lated synthetic M^/Lk values bas ed on the stellar pop- 
ulation models of Maraston (1991) from updated model- 
grids made available by the author 7 and converting from 
their definition of Mk s Q = 3.41). The measurements 
are, therefore, in reasonable agreement with the SSP 
models assuming a Kroupa IMF, given the uncertainties 
in both the data and models. 

Clearly the X-ray determination of the stellar mass 
contribution in these objects can benefit from deeper ob- 
servations, and the systematics involved in the modeling 
of the stellar profile, like the value of the effective radius, 
impact the results considerably. However, the excellent 
agreement between the gravitating mass-to-light ratio at 
the effective radius obtained from X-rays and globular 
cluster kinematics fo r the elliptical galaxy NGC 4649 by 
iBridges et al.l (|2006t ) provides strong support for the reli- 
ability of the stellar mass-to-light ratio determined from 
X-rays in that system. 

Stellar mass has not been detected for NGC 5044, RGH 
80 and NGC 4325, systems which do have Chandra data 
allowing a reasonable sampling of the inner core. Better 
Chandra data would be required in the core of NGC 5044 
to measure possible localized disturbances due to AGN 
activity, which can be a likely source of systematics in the 

http://www-astro.physics.ox.ac.uk/ maras- 
ton/Claudia's_Stellar _Population_Models.html 



mass measurement in the inner 10-20 kpc. Hint of dis- 
turbances (cavities) du e to AGN activity h ave also been 
detected in NGC 4325 (|Russell et al.ll2007t ). There is ev- 
idence from the optical and the Chandra X-ray image 
that the core of RGH 80 may not be completely relaxed. 
Therefore localized departures form hydrostatic equilib- 
rium in the core of these systems are likely explanations 
for the failure to detect stellar mass. 

If we allowed the DM profile to be modified by adi- 
abatic contraction, we obtained substantially smaller 
M*/Lk for our data, which are more discrepant with SSP 
models, casting doubt on the importance of the adiabatic 
contraction process. We ob tained similar results in our 
study on elliptical galaxies (Humph rey et al.l 120061 and 
discus sion therein of ot h er ob servational results). Re- 
cently iGustafsson et alj (|2006l ) proposed that also the 
details of the feedback and its effect upon the concen- 
tration of the baryons is an important ingredient for the 
determination of the final contracted dark matter halo. 

We tested the sensitivity of the measured concentra- 
tions and gravitational masses to the stellar template 
by substituting the effective radius of the De Vaucoleurs 
model derived by 2MASS with the larger value listed in 
RC3. By doing that the stellar mass-to-light ratio can 
increase up to 80% but the concentration and gravita- 
tional mass are not affected as dramatically (see £)9.8p . 
The latter measurements seem therefore robust to the 
assessment of the precise best fit model for the light pro- 
file of BCGs and the precise measuremen t of their size (a 
task which requires part icular care, e.g.. iBernardi et ail 
l2007t ISeigar et alJ I2006L and it is beyond the scope of 
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this paper). 

10.3. Gas fractions 

In Fig. [5] we plot the gas fractions for the objects in 
our sample, calculated within a radius corresponding to 
an over-density of A = 2500 and A = 1250, as a func- 
tion of virial mass. No significant trend of gas fraction 
with mass is present, while there is significant object-to- 
object scatter: the average value of the gas fractions are, 
/ S as,2500 = 0.053 ± 0.012 and / sos ,i 2 50 = 0.069 ± 0.014. 
The mean value of / ga s,2500 obtained for the groups in 
our sample is significantly smaller than that obtained 
fro m the hot, massiv e clusters (T > 5 keV) studied 
by lAh^ri^al] (l2004h . / 3as ,2500 = 0.118 ± 0.016, and 
IVikhlinin etHI (|2006f l. / 9OS , 2 500 = 0.092 ± 0.004. (Note 
we quote the mean and standard deviation, not the gaus- 
sian error- weighted mean and error, because of the pos- 
sibility of non-gaussian contributions to the gas fraction 
distributions, such as intrinsic scatter caused by scale- 
dependent feedback processes.) The fractional error ob- 
tained for our groups, cr/// = 0.2, exceeds t he values of 
0.14 and 0. 04 for, respective l y the clusters of lAllen et al.l 
(|2004l) and IVikhlinin et alj (|2006l) . Therefore, there is 
a clear mass dependence on the gas fraction (mean and 
fractional scatter), not surprising given that the expected 
feedback energy injection by AGN should be more severe 
at the group scale. 

The extrapolation of gas quantities outside the radial 
range of the data is dangerous (see §9.7p . However, for 
most of the groups in our sample, we find that the ex- 
trapolated gas fraction, coupled with the estimate of the 
stellar mass, yield global baryon fractions consistent with 
the universal value; i.e., consistent with notion that X- 
ray b right groups are baryonically closed (| Mathews et al.l 
2005). This result suggests that for the objects in our 
sample for which the slope of the gas density profile is 
not too flat (/3 > 0.5), the extrapolation of the gas frac- 
tion is also fairly reliable. Data at large radii are much 
needed to further explore this issue. 

10.4. Temperature profiles 

The temperature profiles of our groups (Figures 
[31 21 03 an d El) exhibit the sam e behavior char- 
acter i stic of cool core cluste r s (|Markevitch et al.1 
1991 IPeGrandi fe Molendil [200l JPiffaretti et all 120051: 
Vikhlinin et alj l2005t iPratt et al I l2006h : i.e.. the tem- 



perature profile rises outwards from the center, reaches 
a maximum, and then falls at large radius. To exam- 
ine the self-similarity of the profiles, we first rescaled 
them in terms of the virial radii. Then we normalized 
each profile according to the gas-mass-weighted temper- 
ature (T gmw ) computed between 0.1-0.3 r vlr using the 
temperature and gas density models we derived for each 
system. (Gas-mass weighting should be more closely re- 
lated to t he gravitational potential th an emission weight- 
ing, e.g.. iMathiesen fe Evrardll200lD We find that the 
scaled temperature profiles are approximately self-similar 
for r > 0.15r v i r , but there is a large amount of scatter 
at smaller radii (r < 0.1r V j r ). This beh avior is qualita- 
tively similar to that found in clusters by I Vikhlini n et all 
(2006) , though the large scatter in the cores suggest that 
feedback (AGN) processes have had a more dramatic im- 
pact at the group scale; e.g., the striking isothermal core 
of AWM 4. 



Unlike the rise-then-fall temperature profiles observed 
for relaxed, cool core groups and clusters with M > 
10 13 Mq, the three galaxy scale system s (M < 10 13 M ) 
we studied in iHumphrev et al.l (|2006f ) all have temper- 
ature profiles that decrease monotonically with increas- 
ing radius. Hence, the temperature profiles observed in 
our present stu dy provide further suppo rt for the sugges- 
tion we made in (Humphrey et all (|2006l ) that » 10 13 M Q 
represents the mass-scale demarcating the transition be- 
tween (field) galaxies and groups. The dramatic change 
in M/L ratios obser ved at this mass sc ale from optical 
and lensing studies (jParker et al.l I2005L and references 
therein) gives additional evidence that s» 10 13 M Q is a 
special mass scale. 

10.5. Reliability of X-ray mass and concentration 
estimates 

Key sources of systematic errors in the X-ray determi- 
nation of mass and concentration parameters discussed 
in the literature can be listed as: the applicability of hy- 
drostatic equilibrium, a correct interpretation of the tem- 
perature measured by X-ray s atellites (i.e., spectros copic 
versus emission weig hted, e.g. iMazzotta et al.ll2004h and 
the restricted radial range over which the mass is in- 
ferred. 

The results presented in this paper show how mass 
constraints for X-ray bright groups / poor clusters derived 
from good quality XMM and Chandra observations can 
be of the same quality as obtained for hot, massive clus- 
ters. The objects in our sample have been chosen follow- 
ing similar cr iteria for selecting relaxed c lusters for mass 
studi es (e.g.. iPointecouteau et al.l 120051 : [Vikhlini n et al.l 
2006). Indeed, our results show strong support for a sce- 
nario where hydrostatic equilibrium is an excellent ap- 
proximation. The mass profiles inferred from density and 
temperature profiles are in good agreement with the pre- 
dicted quasi-universal NFW profile, and the concentra- 
tion parameters are as expected (i.e., for ACDM) for the 
masses of these objects. The observed trend toward more 
concentrated halos, as expected from numerical simula- 
tions for relaxed halos which have not experienced a re- 
cent major merger, provides further verification our se- 
lection criteria. The fits based on the equation of hydro- 
static equilibrium were able to model complicated tem- 
perature and density profiles (assuming a simple NFW 
DM halo and central stellar component) which would be 
surprising if the gas is significantly out of hydrostatic 
equilibrium. 

The estimate of the real temperature from the pro- 
jected X-ray temperature has been discussed as a source 
of systematic error in X-ray mass estimates. If a spec- 
trum contains several components with different T and 
metallicity the "spectroscopic temperature" derived from 
a single temperatu re fit is biased toward lower temper- 
ature components (jMazzotta et alJl2004t IVikhlinin et all 
2006). But the particular temperature range of 1-3 keV 
explored in this paper is sensitive to the presence of com- 
plex thermal structure, because the Fe-L shell lines com- 
plex is very prominent in the spectrum and able to d is- 
cern different tempera ture components (|Buotel I^OOOd lbl: 
iBohringer et al.l l2002h . If a spectrum contains compo- 
nents with different temperatures, the residuals have 
a characteristic shape, originally noted in ASCA large 
beam spectra as the 'Fe bias", because the Fe abun- 
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Fig. 9. — Observed gas fractions within r2500 (left panel) and within ri250 (right panel) as a function of the virial mass. Errors bars on 
the virial mass have not been shown for clarity of the plot. 



dance derived from such a spectrum is biased low (|Buotd 
l2000d ). Although the Fe abundance is biased low in such 
cases, the inferred average temperature is not biased for 
~ 1 keV systems , as shown in this paper and in Appendix 
of iBuotel (|2000cj) . (Note that such an underestimate in 
the Fe abundance will lead to an overestimate of the gas 
density. But if the bandpass extends down at least to 
pa 0.5 keV, as it does for the A CIS and EPIC data, the 
Fe bias is reduced significantly (Buote 2000a). Further 
reduction in the Fe bias from previous large-beam ASCA 
studies results from using much thinner annuli, especially 
in the core regions exhibiting the steepest temperature 
gradients where the bias would be largest. Consider- 
ing that the Fe-L line emission contributes less than the 
bremsstrahlung over the ACIS and EPIC bandpasses, 
and the gas density varies as the square root of the emis- 
sivity, the reduced Fe Bias leads to a minor overestimate 
of the gas density that is less than the statistical error 
and other systematic errors considered in this paper.) 

What remains to be determined is whether the gas is 
truly single-phase. The use of annuli with the small- 
est possible radial width, allowed by the superb Chandra 
spatial resolution, is crucial to ensure that the multi- 
phase appearance is not simply due to a single phase tem- 
perature gradient. The multiphase appearance of XMM 
spectra extracted in the core of gr oups like NGC 5 33 
(jKaastra et all I2004D and RGH 80 (|Xue et all l200l is 
caused mainly by this reason. This possible error, i.e., 
assuming the gas is single-phase when it is actually mul- 
tiphase, only could affect the innermost regions where 
the temperature gradient is most pronounced. There- 
fore this possible source of systematic in some systems 
can be relevant only for obtaining the most precise mea- 
surement of the stellar mass, but it is unimportant for 
determinations of the halo mass and concentration. 

For the objects in our sample the scale radius is well 
within the radial range covered by the data. Therefore 
a restricted radial range, even for objects for which we 
reach an over-density of 2500, should not be an impor- 
tant source of systematic error in the measurement of 
concentration parameters, as we showed in £19.71 

We also found our results are not sensitive to mild dis- 
turbances related to the presence of a central AGN or not 
fully relaxed dynamical state ( £)9.6p . The derived concen- 



trations and masses are also quite insensitive to errors in 
the shape of the stellar mass profile of the central galaxy, 

(EH). 

It is still very interesting and desirable to obtain data 
at larger radii with offset observations performed by the 
current generation of X-ray observatories, in particular 
XMM, rather than still rely on ROSAT, in particular to 
measure gas masses and gas fract i on (an d to further con- 
strain the total mass). iMol cndi (2004), addressing the 
issue of our ignorance of the outer regions of hot massive 
clusters, pointed out how a combination of reduction of 
particle background to lower levels compared to the cos- 
mic background and the use of a differential measure, 
to improve the knowledge of the actual background in 
the observation, are key to a successful measure of the 
very low sur face brightness regime at large radii. The 
discussion in iMolendl J2004) is focused in particularly in 
measuring the exponential cut-off of the bremsstrahlung 
spectrum of a cluster with T > 3 keV. For the particular 
case of groups/poor clusters, where instead the tempera- 
ture is determined by the Fe-L shell, the clear separation 
of the source component (at least over the radial range 
where the temperature is not declining at values compa- 
rable to the temperature of the soft Galactic background) 
from all the other background components is an effective 
way of making a differential measure, because we know 
both the source and the background. The use of the im- 
proved XMM capabilities in term of collecting area and 
spectral resolution respect to ROSAT will also lift the 
likely important metallicity-density degeneracy (which 
is even more important at large radii, being the group 
emission due to line emission). Furthermore the planned 
future X-ray observatories, like Xeus and Constellation- 
X will have smaller field of view and the mapping of 
the outer regions of nearby systems will be even more 
demanding in terms of observing time. 

11. SUMMARY 

Using Chandra and XMM data we have obtained 
detailed density, temperature and mass profiles of 16 
groups/poor clusters which were selected to be highly 
relaxed systems with the best available data. In sum- 
mary: 

1. The mass profiles were well described by a two com- 
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ponent model: an NFW model for the DM and a 
De Vaucoleurs stellar mass model for 8 objects. For 
objects without adequate sampling in the inner 20 
kpc and for NGC 5044, NGC 4325 and RGH 80 a 
pure NFW model was a good fit of the data. A 
possible explanation for the failure to detect stel- 
lar mass in NGC 5044 and NGC 4325 is localized 
disturbance by AGN activity and for RGH 80 not 
complete relaxation in the core. For objects with 
evidence of stellar mass, the stellar mass-to-light 
ratio in the K band was found to be in approximate 
agreement with simple stellar population synthesis 
models, assuming a Kroupa IMF. 

2. Adopting more complicated models, like introduc- 
ing adiabatic contraction or the recently proposed 
N04 DM profile did not improve the fits. With the 
available data AC produces too low stellar mass- 
to-light ratios and N04 has too high inverse Sersic 
indexes. 

3. The measured c v i r -M v j r relation agrees with the 
predictions of ACDM with a% = 0.9 and f2 m =0.3. 
In particular in the mass range of our group sam- 
ple the expected decrease of c v ; r with M v - lr has been 
detected for the first time. There is a trend, com- 
mon to all X-ray observations, toward more concen- 
trated halos, which can be understood in terms of 
a selection bias, already explored in numerical sim- 
ulations, toward relaxed, earlier forming systems. 

4. The gas fraction measured at an over-density of 
2500 is lower than the one measured for hot mas- 
sive, clusters, and has higher scatter, as expected 
if feedback has played a more severe role at this 
mass scale. However the gas fractions increase 
with radius, and for objects with data extending 
to large radii, these objects are consistent with 
being baryonically closed. However the gas frac- 
tions increase with radius, and for objects with 
data extending to large radii, these objects are 
consistent with being baryonically closed. 

5. When rescaling the radial temperature profiles in 
terms of r v ; r and also the gas-mass-weighted tem- 



perature (evaluated over 0.1-0.3 r v i r ), we find the 
scaled profiles show a fair amount of similarity be- 
yond 0.15 r vir . In the core (r < 0.10 r vir ) the 
scaled profiles have a large amount of scatter, again 
suggesting the important role of feedback in these 
groups. 

We tested the robustness of our results performing 
a careful analysis of possible systematic errors, like 
background subtraction, departures from hydro- 
static equilibrium, deprojection method and found 
none of them to seriously affect our analysis and 
results. 
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APPENDIX 

A: MASSES AND CONCENTRATION PARAMETERS AT DIFFERENT OVER-DENSITIES 



Following lArnaud et all ([2005) and I Vikhlinin et all i|2006h we quote in Table [7] mass, concentration parameters and 
characteristic radii at different over-densities. The virial quantities obtained in our fits have b een rescaled to the ofte n 
used over-density levels of A = 500, 200 and A v ; r listed in Table [1] using the fitting formula of lHu fc Kravtsov! (2003), 
which for c < 20 is accurate to 0.3%. 
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TABLE 7 

Results for the NFW profile at different over-densities 



Group 


Cvir 


f*vir 


M vil 


C500 


»"500 


M 5 oo 


C200 


»*200 


M 200 


(kpc) 




(kpc) 


(10 13 M ) 




(kpc) 


(10 13 M ) 




(kpc) 


(10 13 M Q ) 


NGC 5044 


11.1 ± 0.3 


860 ± 9 


3.74 ± 0.12 


5.7 ± 0.1 


443 ± 4 


2.51 ± 0.07 


8.4 ± 0.2 


653 ± 7 


3.21 ± 0.10 


NGC 1550 


17.0 ± 1.1 


811 ± 13 


3.11 ± 0.15 


9.0 ± 0.6 


428 ± 6 


2.24 ± 0.09 


13.0 ± 0.9 


621 ± 9 


2.73 ± 0.12 


NGC 2563 


9.9 ± 3.4 


762 ± 55 


2.63 ± 0.55 


5.1 ± 1.9 


390 ± 22 


1.72 ± 0.28 


7.5 ± 2.6 


577 ± 39 


2.24 ± 0.43 


Abell 262 


8.9 ± 0.7 


1232 ± 38 


11.00 ± 1.06 


4.5 ± 0.4 


624 ± 15 


7.02 ± 0.52 


6.7 ± 0.5 


930 ± 27 


9.29 ± 0.82 


NGC 533 


16.9 ± 1.2 


727 ± 12 


2.29 ± 0.11 


9.0 ± 0.7 


385 ± 5 


1.65 ± 0.06 


13.0 ± 0.9 


559 ± 8 


2.02 ± 0.09 


MKW4 


12.3 ± 0.8 


1012 ± 18 


6.24 ± 0.34 


6.4 ± 0.5 


527 ± 8 


4.27 ± 0.18 


9.4 ± 0.7 


773 ± 13 


5.40 ± 0.27 


IC 1860 


9.5 ± 0.9 


946 ± 29 


5.07 ± 0.46 


4.9 ± 0.5 


484 ± 11 


3.31 ± 0.23 


7.2 ± 0.7 


718 ± 20 


4.30 ± 0.36 


NGC 5129 


14.6 ± 2.4 


636 ± 29 


1.54 ± 0.22 


7.7 ± 1.3 


335 ± 12 


1.09 ± 0.12 


11.2 ± 1.9 


488 ± 21 


1.35 ± 0.18 


NGC 4325 


11.2 ± 1.2 


833 ± 57 


3.49 ± 0.81 


5.8 ±0.7 


432 ± 25 


2.36 ± 0.45 


8.6 ± 1.0 


635 ± 41 


3.01 ±0.65 


ESO5520200 


7.6 ±0.7 


1288 ±61 


13.03 ± 1.95 


3.9 ±0.4 


650 ± 24 


8.05 ±0.94 


5.8 ±0.6 


976 ± 43 


10.90 ± 1.51 


AWM 4 


8.9 ±0.8 


1384 ± 53 


16.19 ± 1.82 


4.6 ±0.5 


708 ± 23 


10.44 ±0.99 


6.8 ±0.6 


1054 ± 38 


13.75 ± 1.46 


ESO3060170 


8.8 ± 1.0 


1436 ± 141 


18.21 ± 7.58 


4.5 ±0.6 


733 ± 59 


11.67 ±3.72 


6.7 ±0.8 


1093 ± 100 


15.44 ±5.86 


RGH 80 


9.9 ± 1.0 


771 ± 15 


2.81 ± 0.16 


5.1 ± 0.5 


398 ± 5 


1.86 ±0.07 


7.6 ±0.8 


588 ± 10 


2.41 ±0.12 


MS 0116.3-0115 


6.3 ± 2.1 


1272 ± 220 


12.79 ±8.06 


3.1 ± 1.2 


634 ± 85 


7.54 ±3.41 


4.8 ± 1.7 


962 ± 153 


10.55 ±5.92 


Aboil 2717 


6.0 ±0.3 


1432 ±31 


18.39 ± 1.22 


3.0 ±0.2 


710 ± 12 


10.68 ± 0.52 


4.6 ±0.3 


1082 ± 21 


15.10 ±0.89 


RXJ 1159.8+5531 


10.6 ±2.6 


1110 ± 177 


9.08 ± 6.86 


5.6 ± 1.5 


585 ± 73 


6.16 ±3.29 


8.3 ±2.1 


861 ± 127 


7.87 ±5.33 
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B: RELATING X-RAY OBSERVATIONS TO SPHERICAL MODELS OF CORONAL GAS WITH VARIABLE PLASMA 

EMISSIVITY 

We relate the spherical models of density (p„) and temperature (T) of the hot gas in elliptical galaxies, galaxy groups, 
and clusters (^SJ) to the parameters obtained from conventional spectral fitting of X-ray data, such as from Chandra 
and XMM. Since much of the relevant material is scattered throughout the literature spanning over at least 30 years, 
this appendix provides an opportunity to give a self-contained, up-to-date presentation. Our treatment fully accounts 
for radial variations of the plasma emissivity which are often neglected, especially when inferring the gas density. We 
follow the standard procedure where emission from coronal plasma characterized by a single temperature (IT) is fitted 
to the X-ray spectrum extracted from a circular annulus (2D) or, if the data have been deprojected, a spherical shell 
(3D). The density and temperature (and abundance) parameters obtained from such a IT fit are compared to the 
emission- weighted (and projected if 2D) spherical models of p g and T (and abundances). It is assumed the annuli 
are chosen to be sufficiently wide so that the detector point spread function may be neglected. We discuss additional 
weighting by detector responses at the end of this section. 

For any spectral quantity, such as luminosity (erg s~ 4 ), that results from integrating between photon energy E% to 
energy E%, we use the notation, 

l{ae) =L d^ (m) 

to represent the bandpass integration. For a coronal plasma emitting within a volume V, at a temperature Tj with 
chemical abundances Z-p e ,i, Zo y i, Zsi,i, etc. the luminosity is given in xspec by, 

L(AE) t = 4tt [D a (1 + z)] 2 norimA xs (T 4 , Z t ; AE), (B2) 

where we have used the symbol Zi to represent all the abundances, z is the redshift, Da is the angular diameter 
distance (cm), and A xs is the plasma emissivity (erg s _1 cm 3 ) corresponding to the XSPEC implementation of the 
relevant coronal plasma code (e.g., APEC, mekal). To give a precise definition 8 of norm, that accounts for spatial 
variations of the temperature and abundances within Vj, we refer to the volume emissivity (erg s _1 cm~ 3 ) of a coronal 
plasma, 

e(x; AE) — n e (x^nn(x)A(T(x), Z(x)\ AE), (B3) 

where, e.g., T(x) = T(x,y,z), and the plasma emissivity in xspec is related to the conventional definition 9 by, 
A xs = 10 14 A. It is convenient to work in terms of the volume mass density of the gas, p g (x), rather than the volume 
number densities, n e (x) and nu(x), separately. To a very good approximation, for this calculation one may assume a 
fully ionized gas of pure H and He in which case, 

n^ 2 -±^^ n H = i-^, (B4) 

o/i m u op m u 

where p is the mean atomic weight of the gas and m u is the atomic mass unit. (Note that using p = 0.62 corresponding 
to the solar He abundance leads to n e = 1.22njj.) By setting equation |B2| equal to the luminosity obtained by 
integrating equation |B3| over the volume V* and solving for norm; one finds, 



norm,- = 



10- 14 (2 + /*)(4-3a0 



4tt[D a (1 + z)] 2 {bp) 2 m 2 u 



1 



A»(T,,Z i; AE) 



p 2 Jx)A^(T(x),Z(x);AE)d 3 x, (B5) 



Vi 



where norm^ oc J p 2 d 3 x if A xs is constant over the volume Vi. If A xs varies over the volume, then the parameters Ti 
and Zi obtained from fitting a IT coronal plasma model to the spectrum with multiple temperatures and abundances 
will reflect average quantities weighted by the emission profile within the volume and the detector response. We defer 
treatment of the detector response to the end of this section and focus now on the deprojection and projection of 
emission- weighted spherical quantities. 

Deprojection Analysis: If the spectra are deprojected so that the norm^ values refer to spherical shells, then equation 
|B5| can be immediately recast in terms of the weighted square of the gas density, 



{Pi) 



4tt[Da(1 + z)] 2 {hpfmj 
10- 14 (2 + m)(4-30) 

1 3 r- 

A^(T i ,Z i ;AE)(r 3 i -rf_ 1 )J ri _ 1 



(B6) 



v., 

p 2 Jr)A^(T(r),Z(r):AE)r 2 dr, (B7) 



where V, = ^-(rf — rf_i) is the volume of the spherical shell. Note, however, that many deprojection programs like 
PROJCT in xspec and others based on the widely used "onion peeling" deprojection method (jFabian et al.lll98l[ ) give 
normj with respect to the volume, -^(r 2 — rf^) 3 ' 2 , representing the intersection of the 3D shell with the cylinder of 

8 In the XSPEC Users' Manual, norm^ is called K . But the term "norm" is actually displayed for this parameter for the coronal plasma 
models like APEC. 

9 This expression is not given in the XSPEC Users' Manual, but it follows from the definition of norm, (i.e., K) in that manual. 
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the same inner and outer radii (see equation IB12I below) . In this case the integral in equation |B7| must proceed over 
the intersecting volume and therefore depends not only on r. Typically we perform the radial integrations by dividing 
up the shells into 5-10 sub-shells each of constant p g and A xs , where the volume of each sub-shell j is J2k=m(i-i)+i * 

using the definitions associated with equations |Bllj and |BI2| in the projection analysis below. The quantity ' (Pg)t 

is a measure of the average density within the relevant volume (intersecting or total) of the spherical shell. 
The average temperature within the shell i is given by, 

_ T(r)pl(r)A™(T(r), Z(r); AE)r 2 dr 

( r r r ;_ i p2( r )Ax S (T(r),Z(r);AE)r 2 dr ' (B8) 

where, as above, the integration proceeds instead only over the intersecting volume if necessary. We set (T)i — Ti, which 
holds exactly for constant A xs within the shell. (Similarly, for any abundance, such as iron, we set {Z-p c )i = Z-pc.i-) Since 
A xs generally varies monotonically with increasing radius within the shell, this association becomes increasingly more 
accurate as the shell width is allowed to decrease. However, even if calibration uncertainties and other issues associated 
with the detector response can be ignored, systematic errors associated with, e.g., background subtraction, assumption 
of spherical symmetry, Galactic absorption, may bias the inferred IT model fitted to the multi-component spectral 
data over a limited energy range. It is therefore essential to examine the sensitivity of one's analysis to such effects 
as we have done here and previously (|Lewis et al.l l2003; Bu ote fc Le wis 2004; Hu mphrey et al.|[2006h IZappacosta et al.1 
120061) . 

Projection Analysis: This is the primary method employed in this paper. The IT models are fitted directly to the 
spectra that are extracted from concentric circular annuli on the sky. Rather than dividing normj by the emitting 
volume, now we divide it by the area of the annulus, — ir(Ri — i?i_i) 2 , to obtain a line-of-sight, emission-weighted 
projection of pi averaged over Aj, 



proj (p 



47r[D A (l + z)] 2 {bpfmj 
lO" 14 (2 + /i)(4-3/i) 



A 



1 , , 

normi (B9) 



1 2 



A-(T^;AE) {Rj - i? 2 _ x ) J Rt _ 



P RdR [ p 2 (r)A xs (T(r),Z(r);AE)dz, (BIO) 

JRi-i Jlos 



where r = \J R? + z 2 and the projection proceeds within a cylinder denoted by i defined by inner radius, R4—1, outer 
radius, Ri, and height specified by the l.o.s. integration. (Typically we set the l.o.s. integration limits to ss ±1.5 virial 
radii, though our results are very insensitive to this choice for values greater than the virial radius.) For the special case 
of constant A xs within the cylinder i we have that proj (p 2 g )i equals J los Pgdz averaged over A$ (e.g.. lBuote fc Lewisl 
12001 . 

In practice, for fast numerical evaluatio n it is preferable t o approximate the integrations in equation |B10j in terms 
of the contributions from discrete shells (jKriss et al.l Il983h . The spherical volume of the cluster is partitioned into 
a series of N concentric spherical shells such that, ro < r± < T2 < ... < rjv, where the number of shells and 
their widths are chosen to achieve desired computational accuracy. A particular shell, [rfc_i, rfc], is denoted by the 
index fc of the outer radius. We define a corresponding set of N concentric circular annuli on the sky such that, 
Ro — r o < Ri = r i < • • • < Rn = r N, where the origins of the annuli and shells coincide. In analogy to the 3D shells, 
we denote a particular annulus, [Rk-i, Rk], by the index fc of the outer radius. To insure computational accuracy this 
set of reference annuli necessarily over-samples the set of annuli used to extract the X-ray data. Consequently, for 
each annulus i of the data, [iZt— 1, Ri], one defines a mapping between i of the data annuli and k of the reference annuli 
so that the annulus i contains multiple reference annuli. It is useful to consider the case were the reference annuli 
over-sample the data by some integer m so that k = mi; i.e., annulus [Ri-i,Ri] of the data contains all reference 
annuli between R m u-i) and R m i 10 . In this case equation |B10j may be approximated as, 

Proj E E(^) fe V- (Bll) 

V ' ' Z ~ 1J j=m(t-l) + l k=j 

The inner sum X^fclj projects (p 2 A xs ) fc = /9 2 (ffe)A xs (T(ffc), Z(fk)', AE) into reference annulus j, where ffc represents 
an intermediate radius within shell k. The projection is carried out via the matrix, 



47r r 



{rl-R^y -(rl-R]Y + (rl_ l -R)Y -(rl^-R^y , k > j, (B12) 



representing the v olume of shell k that inters ects t he cylinder defined by the width of annulus j and infinite height 
(|Kriss et al.lll983T ). If any terms in equation |B12j have negative arguments they must be set to zero; e.g., if fc = j 
then VjfjJ = Tp(rjfc — r fe-i) 3 ^ 2 ! as noted above in the deprojection analysis. 

10 Since the data do not extend out to the adopted edge of the system, it follows that rjv 2> r m £>, where D is the number of data annuli. 
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The average temperature within annulus i is given by, 

/rp\ SkL RdRf los T(r) P l(r)A™(T(r),Z(r);AE)dz 

(1 )i = 5 Bio 

InL RdR Ls ^(r)A-(T(r), Z(r); AE)dz 

_ S^m(i-l) + l EfcLj ( r Pg AXS )fc Vfc^ . . 

where again we set (T)j = Tj with the same caveats noted above in the deprojection analysis. 

Response Weighting: The existence of radial gradients in the temperature and abundances of galaxy groups and clusters 
imply that the X-ray spectra evaluated over spherical shells or circular annuli of finite width are not IT coronal plasmas. 
In theory, fitting a IT model to such a multi-component plasma yields average values of the spectral parameters norm;, 
Tj, and Zi suitable for comparison with the emission weighted spherical models described previously. Unfortunately, 
biased parameter values can result from such IT fits. Fitting a IT coronal plasma model to a multi-temperature 
spectrum with average tem perature near 1 keV and solar abun dances results in a severe underestimate of the iron 
abundance (" Fe Bias", see iBuote fe Fabianlll998t lBuotdl2000cP . But since the inferred temperature is not biased 
(|Buotd l2000cl) . and the abundance underestimate does not translate to a large overestimate of the gas density, the 
effect on the derived mass profile is minimal. Note that the Fe Bias primarily originates from the differences between 
the strengths of the Fe L shell lines in single and multi-temperature plasmas - not from effects related to the detector 
response other than restricted bandwidth. 

However, since the effective area of the detector responses of ACIS on Chandra and EPIC on XMM (as well as 
detectors on previous X-ray satellites like ROSAT) peak near 1 keV, fitting a IT model to a sp ectrum having a range 
of temperature components above 1 keV will tend to yield a temperature that is biased low (jMazzotta et alJ l2004t 
Vikhlinin 2006^. This effect is most pronounced for the "projection analysis" of galaxy clusters possessing in each 
annulus a wide range of temperatures above ~ 3 keV where the Fe L shell lines are weak. (At lower te mperatures, the 
Fe Bia s ap plies as not e d abo ve.) Approximate methods to account for this effect have been proposed by Maz zotta et alJ 
(|200l and lVikhlinm] (f2006h . 

Our more rigorous approach is a straightforward averaging of the plasma emissivity over the detector response for 
the (projected) region in question. Let the response matrix that determines the probability a photon of energy E 
will be detected in pulse height analyzer (PHA) bin n be denoted by RSPi(n, E), where i denotes the annulus on the 
sky as defined in the "projection analysis" above. This response matrix is usually considered to be the product of 
a "redistribution matrix", RMFi(n, E), and an "auxiliary response file", ARFi(E), the latter of which contains the 
information on the effective area. Then the count rate in PHA bin n is proportional to p 2 J^e A xs (T, Z; £ , )RSP,(n, E), 
where the sum is over all energies in the response matrix. By summing over all PHA bins n corresponding to the 
energy range [E\, E 2 ], and since p 2 does not depend on the convolution, we may account for the detector response in 
our above presentation by replacing A xs with, 

(A»(T, Z- AE)J. = J2H AXS ( T ' Z > £)RSPi(n, E), (B15) 

n E 

where we have not renormalized since the normalization of A xs does not need to be specified for proj (p g )i and (T)j. 

Although applying this response weighting helps to mitigate the temperature bias for hot systems, this procedure is 
still not formally equivalent to that used to obtain 7$ and Zi and norm^ by fitting a IT model to a spectrum containing 
multiple temperature and abundance components. To avoid the Fe Bias and response bias altogether, the projected 
spherical models must be fitted directly to the observed spectra. This means rather than predicting quantities that 
are integrated over the bandpass, one must predict the photon count rate Gi^ n in circular annulus i for each PHA bin 



1CT 14 (2 + m)(4-3m) 



2[D A (l + z)P (5p) 2 m 2 u 

1(T 14 (2 + M )(4-3/x) 
4^[D A (l + z)] 2 (5/i) 2 m2 



f R ' RdR [ dzp 2 g (r)J2^ s (T(r),Z{r);E)RSP l {n,E) 

J Ri _ i J los p 



mi N 

E E 

j=m(i-l) + l k=j 



Vi? (pl) k J2( AXa ( E ))k^i(n,E), (B16) 



such as done, essentially, in a procedure like SMAUG (jPizzolato et al.|[2003[ ). Note that (p 2 )k = p 2 g {rk) and (A xs (i?))fe = 
A xs (T(ffe), Z(fk); E) with fk an intermediate radius in shell k as above. 

It is our experience that even for systems possessing the highest quality data available from Chandra and XMM, the 
magnitudes of statistical errors and the key systematic errors noted above are sufficiently large so that presently there 
is little advantage to increasing the sophistication of the model comparison beyond that expressed by equation |B15| . 
(Note also that we do not use SMAUG in this paper because its current implementation in xspec does not include 
the "parametric mass" fitting approach in <j6l) However, it is expected that data from the next generation of X-ray 
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satellites (Con-X, Xeus) will be of sufficiently high quality to require the direct fitting approach for large numbers of 
systems, provided there is strict control of other systematic errors. 
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